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SUMMARY 


A  method  is  described  for  investigating  the  effects  of  the  dynamic 
response  of  an  autopilot  on  the  stability  of  an  aircraft-autopilot  com¬ 
bination.  The  method  is  based  on  a  study  of  the  constant-damping  curves 
obtained  in  a  plane  defined  by  varying  two  of  the  autopilot  parameters. 
The  dynamics  of  the  autopilot  are  assumed  to  be  describable  by  a  second- 
order  differential  equation.  The  effects  on  the  system  stability  of 
varying  the  gain,  natural  frequency,  and  damping  ratio  of  the  automatic 
damper  are  investigated,  since  these  parameters  determine  the  dynamic 
response  of  the  automatic  damper. 

The  method  is  applied  to  the  analysis  of  the  lateral  motion  of  an 
airplane  equipped  with  a  second-order  automatic  yaw-rate  damper.  For 
any  condition  of  the  airplane,  an  optimum  combination  of  values  of  auto¬ 
pilot  natural  frequency  and  damping  ratio  are  shown  to  exist . for  any 
given  gain  or  required  damping.  A  simple,  analytical  expression  is 
derived  for  obtaining  a  close  approximation  to  these  optimum  points  by 
ignoring  the  effects  of  the  aperiodic  characteristic  modes  of  the  air¬ 
plane.  The  assumption  that  these  aperiodic  modes  may  be  neglected  in 
considering  the  effect  of  the  yaw  damper  on  the  Dutch  roll  oscillation 
is  used  in  all  the  subsequent  analysis.  Expressions  are  derived  for 
the  maximum  damping  obtainable  under  various  conditions.  For  any  given 
natural  frequency  and  damping  ratio  of  the  autopilot,  excessive  auto¬ 
pilot  gain  will  always  cause  the  autopilot  oscillatory  mode  to  become 
unstable. 

Finally,  the  problem  of  designing  an  efficient  yaw  damper  which 
will  improve  the  damping  of  the  Dutch  roll  oscillation  for  various 
flight  conditions  of  an  airplane  is  considered.  A  simple  method  of 
design  is  illustrated  by  applying  it  to  three  flight  conditions  of  an 
airplane.  . 

Calculated  motions,  based  on  the  assumption  of  three  degrees  of 
freedom  for  the  lateral  airplane  motion,  are  presented.  They  agree 
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with  the  results  obtained  from  the  constant-damping-curve  analysis  when 
the  aperiodic  airplane  modes  are  neglected. 


INTRODUCTION 


Recently,  a  great  deal  of  interest  has  been  shown  in  the  use  of 
automatic  stabilization  devices  for  improving  the  damping  of  the  lateral 
(Dutch  roll)  oscillation  of  aircraft  designed  to  travel  at  transonic  and 
supersonic  speeds.  The  analyses  of  such  automatic  stabilization  systems 
may  be  divided  into  two  classes.  In  one  kind  of  analysis  the  effects  of 
various  types  of  autopilots  on  aircraft  stability  are  considered,  and 
usually  the  autopilot  is  assumed  to  be  an  idealized  system  with  no  lags. 
In  the  other  kind  of  analysis  the  effects  of  the  dynamic  response  of  a 
particular  type  of  stabilization  system  on  the  stability  of  the  aircraft- 
autopilot  combination  are  considered. 

Some  investigations  of  the  effects  of  various  types  of  idealized 
autopilots  on  airplane  lateral  stability  are  reported  in  references  1 
to  3-  References  1  and  2  are  theoretical  analyses,  whereas  in  refer¬ 
ence  3  experimental  results  are  compared  with  theoretically  predicted 
effects . 

Analyses  of  the  second  kind,  in  which  the  effects  of  certain  types 
of  dynamic  response  in  a  given  stabilization  system  are  considered,  have 
involved  various  approaches  to  the  problem  of  determining  the  stability 
of  the  complete  system.  The  well-known  frequency-response  analysis  has 
been  used  in  many  papers  (for  example,  in  ref.  4).  Applications  of 
Nyquist's  stability  criterion,  which  was  originally  developed  for  feed¬ 
back  amplifiers,  have  also  been  used.  For  example,  in  reference  5 
Nyquist's  criterion  is  extended  to  systems  with  constant  time  lag  in 
the  feedback  circuit. 

When  the  conditions  which  will  insure  a  given  amount  of  stability, 
rather  than  just  neutral  stability,  are  sought,  the  use  of  constant¬ 
damping  curves  has  been  found  convenient.  The  method  is  described  in 
reference  6,  and  an  example  of  its  application  is  furnished  by  refer¬ 
ence  7-  In  reference  8  a  semigraphical  method  is  developed  for  obtaining 
the  conditions  which  insure  neutral  stability  of  an  airplane-autopilot 
system  when  a  constant  time  lag  in  the  autopilot  is  assumed;  also,  a 
procedure  is  indicated  for  using  the  method  of  constant -damping  curves 
to  determine  the  conditions  which  insure  a  given  amount  of  damping. 

This  method  is  described  in  detail  in  reference  9* 

In  the  present  paper  the  constant-damping- curve  analysis  is  applied 
to  an  airplane-autopilot  system  in  which  the  autopilot  dynamics  are 
represented  by  a  second-order  differential  equation.  Generally,  the 
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actual  frequency  response  of  most  autopilots  (especially  yaw  dampers) 
is  better  approximated  when  a  second-order  differential  equation  is 
used  to  describe  the  dynamic  characteristics  of  the  autopilot  than  when 
the  constant-time-lag  response  is  assumed.  Physically,  this  assumption 
implies  that  the  autopilot  may  be  represented  as  a  damped  oscillatory 
system.  The  present  paper  examines  in  some  detail  the  effects  on  the 
airplane-autopilot  stability  of  variations  in  the  gain,  natural  fre¬ 
quency,  and  damping  ratio  of  the  autopilot. 


A,  B 


SYMBOLS 


coefficients  of  second-order  differential  equation  for  auto¬ 
pilot  dynamics  (see  eq.  (l)) 


wing  span,  ft 


C1  = 


trim  lift  coefficient,  W  cos  y/qS 
rolling-moment  coefficient.  Rolling  moment/qSb 
yawing-moment  coefficient.  Yawing  moment/qSb 
lateral-force  coefficient.  Lateral  force/qS 


8c  7 

8Cy 

C  - 

nP  '  ST 

r  _  Scz 

Zp  8^ 

2V 

cY  =^X 

Yp  8** 

C  - 

np  aEE 

r  = 

-Nrb 

°2v 

8cv 

Cy  =  — X 

1r  ^rb 

°2V 

c 

°2V 

C  -  ^2. 

n5  85 

time-derivative  operator,  d/dt 
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F(D) 

FC(D). 

fn(R)  = 

i  =  fi 
K 

Ko 

% 

KZ 

Kxz 

in 

n 

P,  Q 
Poj  Qo 

p  =  D0 

q 

R 

r  =  Dt 
S 

Tl/2 


characteristic  polynomial  of  airplane-autopilot  system 

characteristic  polynomial  of  airplane  alone 

JL_  dnF(R) 

^  dRn 

gearing  ratio  of  second-order  autopilot,  deg/ deg/ sec 

gearing  ratio  required  for  an  ideal  no-lag  rate  damper  to 
provide  a  given  amount  of  damping,  deg/deg/sec 

nondimens ional  radius  of  gyration  in  roll  about  longitudinal 
stability  axis 

nondimens ional  radius  of  gyration  in  yaw  about  normal  '  ! 
stability  axis 

nondimens ional  product-of-inertia  parameter 

mass  of  airplane,  slugs 

■I. 

an  integer 

coefficients  of  a  quadratic  factor  of  characteristic  equation 
of  airplane-autopilot  system 

coefficients  of  Dutch  roll  quadratic  (see  eq.  (10)) 


dynamic  pressure,  ipV^,  lb/sq  ft 
real  part  of  characteristic  root,  sec"P 


wing  area,  sq  ft 

•  •  '  , 

time  for  amplitude  of  oscillation  to  damp  to  one-half  its 
original  value,  sec 


time,  sec 


6 
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Here  6  is  the  deflection  of  the  control  surface  actuated  by  the  auto¬ 
pilot,  and  X  is  the  component  of  airplane  motion  to  which  the  auto¬ 
pilot  is  sensitive.  In  general,  X  may  be  any  combination  of  the  air¬ 
plane  degrees  of. freedom  or  their  time  derivatives,  or  both.  When  B 
is  positive,  as  it  must  be  for  stable  autopilots,  let  A  =  2£o>0  and 

B  =  a>0^.  Then  equation  (l)  becomes 

D2S  +  2£o)0DS  +  t»026  =  Kid02X  (2) 

Although  A  and  B  are  often  more  convenient  to  use  for  purposes  of 
analysis,  the  equivalent  parameters  £  and  o)Q  are  more  meaningful 

physically.  Most  of  the  results  of  this  paper  will  therefore  be  pre¬ 
sented  in  terms  of  these  parameters.  The  parameter  a>0  is  the  natural 

(undamped)  frequency  of  the  autopilot  system  and  the  parameter  £  is 
the  damping  ratio,  that  is,  the  ratio  of  the  actual  damping  of  the 
system  to  the  critical  damping.  The  parameter  K  which  appears  in 
equations  (l)  and  (2)  is  variously  called  the  amplification  factor,  the 
gain,  or  the  gearing  ratio.  An  autopilot  which  may  be  represented  by 
equation  (l)  or  (2)  is  often  called  a  second-order  autopilot,  since  the 
dynamics  of  the  autopilot  are  represented  by  an  expression  involving 
time  derivatives  of  the  control  deflection  up  to  and  including  second- 
order  derivatives. 

The  coefficients  of  the  characteristic  equation  of  an  airplane 
equipped  with  a  second-order  autopilot  are  functions  of  the  stability 
derivatives  and  mass  characteristics  of  the  airplane  and  of  the  three 
autopilot  parameters  K,  03o,  and  £  (or  K,  A,  and  B) .  These  coef¬ 
ficients  are  often  called  the  stability  coefficients.  If  the  stability 
derivatives  and  mass  characteristics  of  the  airplane  are  known  for  a 
given  flight  condition,  and  if  one  of  the  autopilot  parameters  is 
assigned  some  reasonable  value,  then  the  stability  coefficients  of  the 
airplane-autopilot  system  are  functions  of  the  remaining  two  autopilot 
parameters  only.  Curves  of  constant  damping  of  the  characteristic  modes 
of  the  total  system  may  therefore  be  drawn  in  the  plane  defined  by  con¬ 
sidering  these  two  parameters  as  independent  variables.  Methods  of 
obtaining  such  curves,  with  particular  application  to  oscillatory  modes, 
are  discussed  in  references  6  and  7-  From  an  examination  of  the  constant¬ 
damping  curves  the  ranges  of  values  of  the  two  independent  autopilot 
parameters  which  will  provide  a  given  amount  of  damping  to  the  oscilla¬ 
tory  modes  of  the  airplane-autopilot  system  can  be  determined.  Moreover, 
a  detailed  study  of  these  constant-damping  curves  yields  a  considerable 
amount  of  insight  into  some  of  the  fundamental  properties  of  the  motions 
of  airplanes  equipped  with  second-order  automatic  stabilization  systems. 
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Application  to  Second-Order  Yaw-Rate  Damper 

Equations  of  motion.-  The  method  has  been  applied  to  the  analysis 
of  the  effect  of  a  yaw-rate  damper  on  the  lateral  stability  of  an  air¬ 
plane.  When  the  dimensional  time -derivative  operator  is  used,  the 
lateral  equations  of  motion  are 


-  l  <V  +  -  l  CL^  -  l  % 


2  cnpP  +  2t1bKz''I)2^  -  |  ^  cnrDllf  +  2^bKXZD^  -  r  r-  Cn  = 


2  b  up 


Yl  c  8  V2  c 
b2  Cn55  ^2  Cn 


2ptKxzD^  -  ±  1  CZrDt  +  2pbKx2D2!2(  -  \  1  C2  itf  =  Cz  (5) 


2  b  6P 


In  these  equations  the  assumptions  are  made  that  Cy  =  Cy  = 

and  that  the  only  result  of  the  control  deflection  6  is  a  yawing 
moment.  The  equation  of  motion  of  the  yaw-rate  damper,  written  in 
the  form  of  equation  (2),  is 


u  5  +  2£a>cjD5  +  0)o2S  =  Kb02D^ 


The  complete  characteristic  equation  for  the  equations  of  motion 
(3)  to  (6)  is  a  sixth-degree  equation.  In  fact,  if  Fg(d)  be  the  usual 
fourth-order  characteristic  polynomial  for  the  lateral  degrees  of  free¬ 
dom  described  by  equations  (3)  to  (5)  (see  ref.  2),  then  the  complete 
characteristic  equation  for  the  airplane-autopilot  combination  is 

F(D)  -  (jl-  +  +  c%2)f0(D)  -  j£|  CngKo>02  "jWDB  -  j  Cy  J ^nbKx2D2  - 


Hv)- 


y3 


(7) 
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The  characteristic  modes  obtained  from  this  equation  are  usually  two 
oscillations  and  two  aperiodic  modes.  In  most  cases  the  two  aperiodic 
modes  and  one  of  the  oscillations  may  be  associated  with  the  lateral 
motion  of  the  airplane,  whereas  the  second  oscillation  may  be  associ¬ 
ated  with  the  motion  of  the  autopilot  system. 

Choice  of  variable  autopilot  parameters  for  constant -damping 
curves.-  If  all  the  parameters  in  equations  (3)  to  (5)  are  known  for 
an  airplane  in  a  given  flight  condition,  the  effects  of  K,  (Dq,  and  £ 
on  the  lateral  oscillatory  stability  of  the  combined  system  Can  be 
determined  by  varying  these  parameters  in  equation  (7),  as  is  shown  in 
appendix  A.  One  of  these  parameters  must  be  fixed  if  the  stability 
boundaries  are  to  be  plotted  in  a  plane  as  described  in  appendix  A,  and 
therefore  the  relative  convenience  of  fixing  each  of  the  three  param¬ 
eters  should  be  considered.  Since  the  dynamic  characteristics  of  the 
autopilot  may  be  expressed  in  terms  of  its  frequency  response,  a  study 
of  the  effects  of  the  three  autopilot  parameters  on  this  frequency 
response  should  provide  some  insight  into  the  question  of  which  two 
parameters  should  be  varied  simultaneously., 

The  transfer  function  of  the  stabilizing  autopilot  may  be  obtained 
from  equation  (6): 

O 

r-  ,  -1  KOU 

[s/D+IpCD)  -  — - 2— - -  (8) 

D2  +  2£0>oD  +  CD02 


Therefore,  the  expression  for  the  autopilot  frequency  response  is 


(&/Dtj  (in>)  = 


K 


Rp(co) 


10p(cn) 


(9) 


1  _  +  2i5  £- 

o>  2  “o 


The  frequency  response  of  the  autopilot  is  obtained  from  equa¬ 
tion  (9)  by  plotting  the  amplitude  Rp(o>)  and  phase  angle  0p(ao)  of 

the  complex  number  [s/D\j/]p(  io>)  against  cd.  The  phase  angle  is  inde¬ 
pendent  of  K,  whereas  the  magnitude  at  any  value  of  CD  is  proportional 
to  K.  Thus  the  phase-angle  curve  and  the  shape  of  the  magnitude  curve 
both  depend  only  on  £  and  <jdq,  and  the  gearing  ratio  simply  acts  as  an 
amplification  factor  on  the  magnitude  curve.  From  this  point  of  view, 
it  would  seem  desirable  to  select  reasonable  values'  of  K  and  allow  £ 
and  o>0  to  be  the  variable  parameters  in  equation  (7)*  A  study  of  the 
stability  boundaries  in  the  -plane  would  then  show  the  effect  of 
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varying  the  shape  of  the  autopilot  frequency  response  on 'the  stability 
of  the  total  system,  whereas  the  variation  of  K  would  show  the  effect 
of  varying  the  amplification  factor.  The  stability  boundaries  were 
therefore  first  calculated  in  the  £cjDo-plane,  and  some  of  the  effects 
of  varying  K  on  these  curves  were  investigated. 

Since  the  gearing  ratio  K  is  in  many  respects  the  most  important 
parameter  determining  the  stability  of  the  total  system,  stability  bound¬ 
aries  were  also  obtained  in  the  Kh)0-plane  for  fixed  values  of  £ .  This 
method  of  plotting  the  stability  boundaries  has  two  distinct  advantages. 
First,  the  two  most  important  parameters,  K  and  cn0,  are  allowed  to 
vary.  Second,  the  fixed  parameter  £  is  known  to  be  between  0  and  1 
for  stable  oscillatory  autopilot  systems;  therefore,  the  effect  of 
varying  £  on  the  stability  boundaries  in  the  Kh>0-plane  can  be  deter¬ 
mined  fairly  easily  by  choosing  several  values  of  £  which  will  span 
this  range. 

Constant-damping  curves  with  gearing  ratio  fixed;  three  degrees  of 
freedom.-  In  order  to  investigate  the  effect  of  the  second-order  yaw 
damper  on  the  stability  of  the  system,  values  of  the  airplane  parameters 
were  inserted  into  equations  (3.)  to  (5)  to  correspond  to  a  cruising 
flight  condition  for  the  airplane  described  in. table  I.  The  Dutch  roll 
oscillation  for  this  flight  condition  has  a  period  of  1.30  seconds  and 
T1/2  of  2.6o  seconds.  1 

In  order  to  investigate  the  type  of  constant-damping  curve  which 
appears  in  the  £<n0-plane  (or  AB-plane)  for  constant  K,  a  value  of  K 
was  chosen  which  would  give  good  damping  if  used  in  a  perfect  propor¬ 
tional  yaw  damper  that  has  no  inertia  or  damping.  The  value  K  =  0.086 
degree  of  rudder  deflection  per  degree  per  second  of  yawing  velocity 
was  chosen,  which  would  make  T^y^  Dutch  roll  oscillation  equal 

to  0.75  second.  The  constant-damping  curves  in  the  AB-plane  were  then 
drawn  for  this  value  of  K  by  using  the  equations  given  in  appendix  A, 
and  are  shown  in  figures  1  to  5«  The  frequencies  of  the  modes  on  a  . 
given  curve  vary  from  zero  to  infinity,  and  typical  values  are  shown  on 
the  individual  curve  s . 

As  pointed  out  in  appendix  A,  it  is  simpler  to  obtain  the  curves 
first  in  terms  of  A  and  B  and  then  substitute  the  more  significant 
parameters  £  and  co0.  However,  for  a  preliminary  investigation  of 
the  general  types  of  curves  and  the  manner  in  which  they  change  as  R 
and  cjd  (the  damping  and  frequency  parameters)  are  varied,  continued 
use  of  the  AB-plane  is  more  convenient  because  negative  values  of  B 
may  be  considered,  whereas  the  parameters  £  and  0>o  have  an  obvious 
physical  significance  only  for  positive  values  of  B.  For  this  reason,- 
the  typical  curves  of  figures  1  to  5  are  drawn  in  the  AB-plane.  These 
figures  are  presented  and  discussed  primarily  to  familiarize  the  reader 
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with  the  general  types  of  curves  which  occur,  so  that  the  later  com¬ 
parison  of  the  curves  in  the  £ti>0_Plane  will  not  be  too  confusing.  The 
practically  significant  portions  of  the  curves  are  essentially  the  same 
in  either  plane. 

Figure  1  presents  the  zero-damping  curve >  which  is  continuous  and 
crosses  itself  to  form  a  loop.  The  loop  has  no  practical  significance, 
since  it  merely  defines  the  region  of  autopilot  parameters  for  which 
there  are  two  unstable  oscillations.  The  hatch  marks  indicate  the 
boundary  of  the  region  in  which  all  modes  are  stable.  In  all  cases  the 
B  =  0  axis  is  the  boundary  at  which  an  aperiodic  mode  becomes  unstable. 

The  curve  of  greatest  practical  importance  is  that  defined  by  the 
damping  of  the  airplane  without  the  autopilot;  that  is,  the  curve  for 
Tp/2  =2.6o  seconds.  This  curve,  shown  in  figure  2,  is  discontinuous. 

When  a)  approaches  the  airplane  frequency  (the  frequency  of  the  Dutch 
roll  oscillation  of  the  airplane  alone) ,  both  A  and  B  become  infi¬ 
nite  in  magnitude.  These  infinities  are  caused  by  the  vanishing  of  a 
factor  in  the  denominator  of  the  expressions  for  A  and  B.  At  the 
airplane  frequency  a  new  branch  of  the  constant-damping  curve  is  started. 
The  region  in  the  AB-plane  in  which  there  are  no  oscillatory  modes  that 
have  less  damping  than  the  airplane  alone .(indicated  by  hatch  marks)  is 
bounded  by  this  new  branch  of  the  curve.  The  dashed  line  shows  the 
value  B  =  Qq  =  23.84,  where  Q0  is  the  constant  coefficient  of  the 

Dutch  roll  quadratic,  which  is  written  as  D^  +  P0D  +  Q0.  The  signifi¬ 
cance  of  this  value  will  be  discussed  subsequently. 

Figure  3  shows  a  typical  curve  for  a  damping  somewhat  greater  than 
that  of  the  airplane  alone.  It  has  two  points  of  discontinuity  and 
three  branches.  Thus,  the  airplane -damping  curve  shown  in  figure  2  is 
a  critical  curve,  separating  the  continuous  curves  for  less  damping 
from  the  doubly  discontinuous  curves  for  greater  damping.  The  two 
critical  values  of  cd  are  very  close,  so  that  the  branch  of  the  curve 
in  the  negative-B  region  corresponds  to  a- very  small  range  of  co. 

This  branch  of  the  curve  is  of  academic  interest  only.  The  part  of  the 
curve  in  figure  3  which  forms  the  significant  boundary  is  again  the 
curve  in  the  upper  right  quadrant  of  the  plane,  defining  a  region  pf 
better  damping  indicated  by  hatch  marks. 

The  curve  in  figure  4  is  for  ^1/2  =  0-75  second,  the  damping 

which  the  airplane -autopilot  oscillation  would  have  if  the  autopilot 
had  no  lags.  This  curve  has  another  branch  at  large  negative  values 
of  B  and  negative  A,  but  this  branch  is  not  shown  because,  as  has 
been  mentioned,  it  is  of  academic  importance  only. .  (it  should  be 
noted  that  the  significant  portions  of  all  these  discontinuous  curves 
start  at  values  of  03  near  the  airplane  frequency.)  The  curve  for 
T1/2  =0.75  second  also  represents  a  critical  damping  value,  since 
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for  larger  values  of  damping  the  curves  again  become  continuous  and 
form  loops,  as  shown  in  figure  5  for  T =  O.60  second.  In  the 

following  discussion  the  damping  corresponding  to  Tp/g  =  0*75  second 

and  the  associated  constant -damping  curve  will  be  called  the  critical 
damping  and  critical  curve.  ■  1 

These  curves  will  be  discussed  in  greater  detail.  For  the  present 
it  will  suffice  to  note  that,  for  any  gearing  ratio,  greater  damping 
can  be  provided  by  a  second-order  autopilot  than  could  be  provided  by 
a  perfect  proportional  autopilot  with  the  same  gearing  ratio,  provided 
that  the  proper  values  of  £  and  u>0  are  chosen  from  a  closed-loop 
region  of  the  type  shown  in  figure  5* 

Equivalent-oscillator  concept.-  Appendix  A  shows  that,  in  deter¬ 
mining  the  effect  of  the  yaw  damper  on  the  stability  of  the  Dutch' roll 
oscillation,  the  airplane  is  represented  by  an  equivalent  oscillator 
whose  period  and  damping  are  those  of  the  Dutch  roll  oscillation. 
Explicitly,  the  assumption  is  made  that  the  three-degree-of-freedom 
equations  of  motion  of  the  airplane,  given  by  equations  (3)  to  (5), 
may  be  replaced  by  the  single  equation  of  motion 

2^(kz2  -  +  PoD+  +  V)  -  ^  cn6B  =  ^  cn  <10) 

The  physical  interpretation  of  this  assumption  is  that  the  effect  of 
the  aperiodic  modes  on  the  Dutch  roll  stability  is  small.  The  quadratic 

p 

equation  D  +  PQD  +  Q0  =  0  yields  the  complex  characteristic  root 
corresponding  to  the  Dutch  roll  oscillatory  mode. 

The  analysis  is  greatly  simplified  by  this  equivalent-oscillator 
concept,  and  it  will  be  seen  that  this  simplified  analysis  gives  adequate 
results.  Except  where  otherwise  specified,  the  equivalent  oscillator  is 
assumed  to  represent  the  airplane  in  all  the  subsequent  discussions.  The 
curves  in  figures  1  to  5,  obtained  from  the  three-degree-of-freedom 
analysis,  were  drawn  primarily  as  a  check  of  the  accuracy  of  the 
equivalent-oscillator  approximation. 

In  figure  6  the  constant-damping  curves  calculated- by  the  equivalent- 
oscillator  analysis  for  representative  values  of  damping  are  presented  in 
the  £u>0-plane.  In  interpreting  these  curves  it  is  essential  to  remember 
that  two  oscillations  are  present  for  most  of  the  significant  points  in 
the  £o)0-plane.  One  of  these  oscillations  may  generally  be  associated 

with  the  airplane  and  the  other  with  the  autopilot.  Therefore,  at  every 
point  in  the  most  significant  regions  of  the  £o>0-plane  two  damping  curves 

must  cross  each  other. 
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All  the  curves  previously  obtained  from  the  three-degree-of-freedom 
analysis  are  presented  in  figure  6,  and  in  addition  the  boundary  of 
equal  roots  and  the  curve  for  Tp/2  =  O.38  second  are  shown.  Each  of 
the  constant- damping  curves  starts  on  the  boundary  of  equal  roots  at 
CD  =  0.  Before  the  general  discussion  of  figure  6  is  presented,  the 
significance  of  the  curve  for  T2/2  =  0*38  second  will  be  explained. 

,  Maximum  damping  with  fixed  gearing  ratio.-  As  mentioned  previously, 
the  constant- damping  curves  for  higher  damping  than  the  critical  damping 
are  continuous  and  form  loops  defining  the  regions  which  insure  greater 
damping.  Since  the  loops  become  smaller  for  larger  values  of  damping, 
it  is  reasonable  to  assume  that  for  some  value  of  damping  the  loop  will 
become  vanishingly  small  -  only  a  cusp  in  the  curve  (see  ref.  9)*  This 
curve  then  corresponds  to  the  maximum  damping  which  can  be  obtained  for 
the  Dutch  roll  oscillation  with  a  second-order  yaw-rate  damper  for  the 
given  gearing  ratio.  Thus  the  cusp  point  may  be  considered  the  optimum 
point  in  the  f;cn0-plane  for  a  given  gearing  ratio,  and  this  point  deter¬ 
mines  the  shape  of  the  autopilot  frequency  response  which  will  give  the  , 
highest  damping  to- the  Dutch  roll  oscillation  for  a  given  amplification  , 
factor  K. 


These  optimum  points  may  be  obtained  by  a  rather  simple  algebraic 
analysis  because  they  correspond  to  a  double  oscillatory  root  of  the 
characteristic  equation,  which  is  a  quartic  equation  when  the  equivalent- 
oscillator  concept  is  used.  The  derivation  of  the  optimum-point  char¬ 
acteristics  is  presented  in  appendix  B.  The  damping  and  frequency  of 
the  double  oscillatory  mode  are  assumed  to  be  given  by  the  characteristic 
equation  - 

(d2  +'  PD  +  q)2  =  0  (11) 

As  shown  in  appendix  B,  the  value  of  Q  may  be  obtained  by  solving  the 
quadratic  equation,  for  a  given  gearing  ratio, 


C-iK  \  2 

Po  ±  2^/Q 


-  2QqQ  +  Q02  =  0 


(12) 


For  positive  values  of-  K  the  larger  real  root  in  equation  (12)  is  used. 
With  this  value  for  Q,  P  may  be  obtained  from  the  expression 


(p0  +  ClK)<?  -  PoQ0g 
2«o(«  '  Qo) 


(13) 
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These  values  of  P  and  Q,  when  used  in  equation  (ll),  give  the  charac¬ 
teristic  root  (and  therefore  the  period  and  damping)  corresponding  to  the 
best-damped  Dutch  roll  motion  obtainable  with  the  given  gearing  ratio. 

In  order  to  find  the  autopilot  parameters  A  and  B  (or  £  and  cno) 
which  yield  this  maximum  damping,  these  values  of  P  and  Q  are  sub¬ 
stituted  into  the  following  expression's  for  A  and  B: 


(14) 


(15) 


The  maximum  damping  for  the  airplane  under  discussion  for 
K  =  0.086  deg/deg/sec  was  obtained  in  this  manner  and  was  found  to 
.correspond  to  T^/2  =0.38  second.  As  can  be  seen  from  figure  6>  this 
curve  does  have  a  cusp  at  the  optimum  point  in  the  £<r>0- plane . 


Since  the  high-damping  curves  in  the  three-degree-of-freedom 
analysis  also  form  loops,  the  same  type  of  analysis  can /be  used  to 
obtain  the  cusp  point  in  the  {[(ng-plane.  The  calculation  of  the  maximum 
damping  would  be  much  more  complicated,  however,  for  the  three-degree- 
of-freedom  case.  To  show  that  the  maximum  damping  as  calculated  by  the 
equivalent-oscillator  analysis  is  an  adequate  approximation  to  the  maxi¬ 
mum  damping  for  the  complete  airplane,  the  curve  for  T^/g  =0.38  second 

was  drawn  for  the  three-degree-of-freedom  case  also. 


The  comparison  of  the  curves  for  T^2  =  0.38  second  for  the  two 

cases  i.s  shown  in  figure  7.  In  this  figure  the  significant  portions  of 
other  typical  damping  curves  are  also  shown.  For  practical  purposes, 
the  equivalent-oscillator  analysis  is  an  adequate  approximation  to  the 
three-degree-of-freedom  analysis  in  determining  the  required  constant¬ 
damping  curves,  including  the  maximum- damping  point.  It  might  also  be 
noted  that  the  critical  damping  is  not  exactly  the  same  in  the  two  cases 
(T1/2  '=  0.75  second  for  the  three-degree-of-freedom  case  and 

T1/2  =  0.73  second  for  the  equivalent  one-degree-of -freedom  case) . 

Discussion  and  interpretation  of  damping  boundaries  in  £u>Q- plane . - 
Figure  7  shows  that  there  is  little  change  in  the  stability  boundaries 
when  the  lateral  motion  of  the  airplane  is  represented  by  the  equivalent 
oscillator.  This  result  means  that  the  real  characteristic  modes  of  the 
lateral  motion  can  be  neglected  in  calculating  the  effect  of  the  auto¬ 
matic  yaw  damper  on  the  Dutch  roll  oscillation.  The  only  fundamental 
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difference  occurs  in  the  airplane-damping  curve,  Tq/2  =  2.60  seconds. 

For  the  equivalent  oscillator,  this  curve  is  discontinuous  at  a  value  of 
cd0  equal  to  the  airplane  frequency,  4.8  radians  per  second.  The  slope 
of  the  curve  becomes  infinite  at  this  value  of  0>o,'and  the  portion  of 

the  curve  bounding  the  higher-damping  region  begins  here  at  £  =  00 .  In 
the  three-degree-of-freedom  analysis,  the  curve  for  Tq/p  =  2.60  seconds 

has  a  slope  which  becomes  very  large  at  cjo0  =  4.8  radians  per  second 
but  remains  finite.  Thus  the  difference  in  the  curves  is  negligible 
for  practical  purposes.  Since  the  significant  portion  of  the  airplane¬ 
damping  curve  begins  at  o)0  =  4.8  radians  per  second  for  any  positive 
gearing  ratio,  the  damping  of  the  equivalent  oscillator  cannot  be 
improved  with  a  positively  geared  yaw-rate  damper  which  has  a  natural 
frequency  less  than  the  frequency  of  the  oscillator  itself.  Figure  7 
shows  that  for  practical  purposes  the  same  statement  can  be  made  for 
the  actual  airplane. 

Certain  general  effects  of  varying  the  autopilot  damping  and  natural 
frequency  on  the  stability  of  the  system  may  be  observed  in  figures  6 
and  7.  The  area  of  interest  is  the  roughly  rectangular  region  indicated 
by  the  hatch  marks  in  figure  2  in  the  AB-plane,  in  which  both  oscillatory 
modes  have  more  damping  than  the  airplane  alone .  The  value  of  B  indi¬ 
cated  by  the  dashed  line  in  figure  2  corresponds  to  the  airplane  fre¬ 
quency  (u>0  =  4.8  radians  per  second)  .  If  £  is  fixed  at  any  positive 
value  and  cn0  is  increased,  figure  6  indicates  that  the  system  damping 
increases  to  a  maximum  value  at  some  value  of  cjd0  and  then  -drops  off, 
approaching  the  critical  damping  as  tn0  approaches  infinity.  Similarly, 
at  any  fixed  natural  frequency  greater  than  the  Dutch  roll  frequency,  if 
the  value  of  £  is  increased  from  zero  the  damping  reaches  a  maximum  at 
some  value  of  £  and  then  drops  off,  approaching  the  airplane  damping 
as  £  approaches  infinity.  Thus,  for  fixed  ,  £  there  is  an  optimum  co0 , 
and  for  fixed  o>0  an  optimum  £.  The  best  of  all  these  points  is  the 
max imum- damping  point,  which  is  obtained  by  the  simple  calculation  pre¬ 
viously  described. 

Importance  of  oscillation  frequencies  in  interpreting  constant¬ 
damping  curves.-  From  stability  considerations  alone,  the  regions  defined 
by  the  hatch-marked  portions  of  the  curves  shown  in  figures  1  to  5  deter¬ 
mine  the  values  of  autopilot  parameters  which  guarantee  at  least  the 
indicated  amount  of  damping.  However,  points  may  be  chosen  outside  a 
given  region  which  still  seem  to  give  an  airplane  motion  that  is  as  well- 
damped  as  that  for  points  in  the  region.  For  example,  figure  8  shows 
the  motions  obtained  with  autopilots  having  the  characteristics  defined 
by  points  1,  2,  and  3  in  figure  6.  All  motions  shown  in  this  paper  were 
obtained  from  a  Reeves  Electronic  Analog  Computer,  by  use  of  the  three- 
degree-of-freedom  equations  of  motion.  The  fact  that  these  motions  check 
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the  predictions  made  from  the  equivalent-oscillator  stability  boundaries 
confirms  the  adequacy  of  the  equivalent-oscillator  analysis. 

Points  2  and  3  are  on  the  good  portion  (loop)  of  the  curve  for 
lp/2  =  O.60  second,  but  point  1  is  on  the  intersection  of  the  curves 

for  T1/2  =0.60  second  and  I1/2  =  2.60  seconds  and  is  outside  the 

good  region  defined  by  the  loop.  However,  the  actual  airplane  motions, 
represented  by  the  sideslip  and  roll  motions  in  figure  8,  are  very 
similar  for  all  three  autopilots,  and  the  lightly  damped  mode  is  impor¬ 
tant  only  in  the  rudder  motion  for  the  case  corresponding  to  point  1 

(fig.  8(a)).  Actually,  the  effect  of  the  lightly  damped  mode  can  be 

seen  in  the  sideslip  motion  of  figure  8(a),  but  it  is  almost  negligible. 
Although  point  1  in  figure  6  corresponds  to  T-jy2  =  2.60  seconds,  as  far 

as  the  airplane  motion  is  concerned  this  autopilot  would  seem  to  give  as 
good  damping  as  autopilots  whose  characteristics  fall  in  the  loop  of 
T]_^2  =  0.60  second. 

In  order  to  understand  why  the  lightly  damped  mode  corresponding 
to  point  1  in  figure  6  has  practically  no  effect  on  the  airplane  motion, 
the  frequencies  of  the  modes  which  are  predicted  at  this  point  must  be 
considered.  Since  'for  the  most  significant  points  in  the  £co0-plane  there 
must  be  two  characteristic  oscillations,  each  of  these  points  is  actually 
a  crossing  point  of  two  constant-damping  curves,  as  can  be  seen  for 
point  1  in  figure  6.  The  general  trends  of  the  frequencies  along  the 
damping  curves  are  shown  in  figures  1  to  5.  Along  the  final  portion  of 
each  curve  (that  is,  the  portion  which  approaches  the  (7=0  axis),  the 
frequencies  correspond  to  the  autopilot  frequency.  At  point  1  in  fig¬ 
ure  6  the  mode  on  the  curve  for  T1/2  =  2.60  seconds  has  a  frequency 
co  w  10  radians  per  second,  a  value  which  corresponds  to  the  autopilot 
frequency  (as  can  be  seen  from  fig.  2),  whereas  the  better-damped  mode 
is  the  airplane  mode  with  co  5  radians  per  second  (see  fig.  5). 

Since  the  frequency  of  the  autopilot  mode  is  approximately  twice  that 
of  the  airplane  mode,  the  effect  of  the  corresponding  lightly  damped 
rudder  oscillation  on  the  airplane  motion  is  small,  for  the  airplane 
cannot  follow  such  rapid  oscillations. 


Consideration  of  the  frequencies  which  occur  at  points  along  the 
damping  curves  is  thus  seen  to  be  important  in  attempting  to  predict 
the  type  of  motion  which  would  be  obtained  with  autopilots  whose  char¬ 
acteristics  are  determined  by  these  points.  This  frequency  effect  is 
brought  out  even  more  strongly  by  figure  9,  which  shows  the  motions 
corresponding  to  points  4  and  5  in  figure  6.  In  this  case  both  points 
are  on  the  zero-damping  curve.  At  point  5  both  modes  have  approximately 
the  same  frequency.  However,  at  point  4  the  autopilot  frequency  is  the 
neutrally  damped  one,  with  co  **  13. 5  radians  per  second,  and  the  airplane 
mode  lies  on  one  of  the  well-damped  loop  curves,  which  crosses  the  curve 
for  T^/2  =  00  at  point  4.  The  motions  shown  give  the  results  predicted 
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by  the  analysis.  As  shown ’in  figure  9(a),  corresponding  to  point  4,  the 
neutrally  damped  high-frequency  rudder  motion  has  little  effect  on  the 
well-damped  airplane  motion.  On  the  other  hand,  for  point  5  the  effect 
of  the  neutrally  damped  mode,  which  is  close  to  the  natural  airplane 
frequency,  is  the  dominant  neutrally  damped  airplane  oscillation  shown 
in  figure  9(b).  Autopilots  which  cause  very  poorly  damped  control 
motions,  however,  would  be  unsatisfactory  from  a  practical  point  of 
view  even  if  their  effect  on  the  airplane  damping  were  satisfactory. 

Effect  of  varying  K  on  curves  in  -plane .  -  The  effect  of  gearing 
ratio  may  be  obtained  by  considering  the  effect  of  varying  K  on  the 
curves  in  the  £cDQ-plane.  By  using  the  equivalent-oscillator  analysis, 

the  critical  damping  is  found  to  be  Rcr  =  "^(^0  +  ^1^)‘  Since  the  air¬ 
plane  damping  is  -7po>  the  critical  damping  becomes  the  airplane  damping 

as  K  vanishes.  This  result  is,  of  course >  necessary,  since  zero 
gearing  ratio  implies  no  autopilot.  As  K  is  made  smaller,  the  whole 
set  of  loop  curves  in  figure  6  tends  to  move  to  the  left,  since  the 
critical  curve  approaches  the  airplane  damping  curve.  Conversely,  as  K 
increases,  these  curves  move  to  the  right.  Also,  the  loop  corresponding 
to  any  given  damping  larger  than  the  critical  damping  must  expand  as^  K 
increases,  since  the  given  damping  comes  closer  to  the  critical  damping. 
Therefore  the  given  loop  approaches  the  infinite  loop  asymptotic  to  the 
critical-damping  curve.  Physically,  this  result  simply  means  that  as 
the  gearing  ratio  is  increased  there  is  a  larger  range  of  values  of  £ 
and  u)0  for  which  a  given  damping  larger  than  the  critical  damping  may 

be  obtained. 

A  clearer  idea  of  the  way  in  which  the  set  of  loop  curves  moves  in 
the  £cn0 -plane  as  K  varies  is  obtained  by  investigating  the  variation 
of  the  maximum-damping  point  as  K  varies.  The  position  of  the  maximum¬ 
damping  point  is  itself  important,  because  it  is  the  optimum  combination 
of  £  and  cd0  for  any  value  of  Kj  but,  also,  since  this  point  is  a 
kernel  which  is  surrounded  by  all  the  loops,  the  motion  of  this  point 
gives  a  clearer  idea  of  the  motion  of  any  loop  as  K  varies.  The 
desired  variation  may  be  easily  obtained  by  inserting  values  of  K  into 
equations  (12)  and  (13)>  and  using  the  resulting  values  of  P  and  Q 
in  equations  (l4)  and  (15)-  However,  a  simpler  and  clearer  method  is 
shown  to  be  possible  in  appendix  B,  wherein  a  value  of  damping  is  assumed 
and  solutions  are  found  for  the  values  of  K,  A,  B,  and  Q  which  will 
make  this  the  maximum  damping.  This  procedure  clearly  gives  the  smallest 
magnitude  of'  K  with  which  the  desired  damping  may  be  obtained  with  the 
second-order  automatic  yaw  damper,  and  the  associated  values  of  A  and 
B  ' then  may  be  considered  as  "optimum"  values  for  the  given  airplane  and 
desired  damping.  In, this  method  equation  (13)  is  replaced  by 
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where  R  and  T^/2  correspond  to  the  desired  damping.  Equation  (l4) 
then  gives  A  immediately,  and  Q  may  be  obtained  from  the  expression 

Q  =  %  ±  \fo(p  ~  po)  (17) 

Equation  ( 17 )  will  give  two  values  of  Q,  one  greater  than  Qq 
and  the  other  less  than  Q0.  ,  By  using  these  values  in  equation  (15), 
the  corresponding  values  of  B  are  obtained,  and  from  equation  ( 13 ) : 

1  , 

..  2PQ  -  PnB  -  AQn  . 


The  smaller  value  of  Q  obtained  from  equation  (17)  gives  a  value 
of  B  corresponding  to  a  value  of  lower  than  the  airplane  fre¬ 

quency  and  results  in  a  negative  K.  The  larger  value  of  Q  gives  a 
value  of  (n0  higher  than  the  airplane  frequency  and  positive  gearing. 
Figures  10  and  11  show  the  variation  of  the  optimum  points.  Figure  10 
shows  the  curve  on  which  the  points  lie  in  the  £o>0-plane.  Figure  11 
shows  the  gearing  necessary  to  obtain  any  given  damping  as  the  maximum 
damping  (when  the  autopilot  characteristics  are  the  optimum  ones  for 
that  gearing) . 

The  point  corresponding  to  K  =  0  on  each  curve  is  significant 
only  as  a  limiting  point,  since  zero  gearing  implies  no  autopilot. 

Since  the  double  oscillatory  mode  corresponding  to  the  maximum  damping 
must  approach  the  Dutch  roll  mode  as  K  approaches  zero,  the  optimum 
point  in  the  -plane  approaches  the  values  corresponding  to  the  Dutch 
roll  mode.  Increased  damping  of  the  system  can  be  obtained ’with  either 
positive  or  negative  gearing.  The  positive-K  branch  of  the  locus  of 
optimum  points  given  in  figure  10  lies  in  the  range  of  values  of  o>0 
higher  than  the  airplane  frequency.  Thus,  as  mentioned  previously, 
second-order  yaw  dampers  with  positive  gearing  must  have  a  value  of  cd0 
higher  than  the  airplane  frequency  in  order  to  improve  the  damping. 

Also,  the  variation  of  the  optimum  points  on  the  positive-K  branch  of 
figure  10  shows  that,  if  increased  damping  from  the  autopilot  is  sought 
by  increasing  the  gearing  ratio,  the  natural  frequency  and  damping  ratio 
of  the  autdpilot  should  generally  be  increased  simultaneously.  This 
fact  can  be  of  considerable  practical  importance,  as  will  be  brought 
out  more  clearly  in  subsequent  discussion. 

Comparison  of  advantages  of  positive  and  negative  gearing.-  If 
negative  gearings  are  used,  <n0  values  lower  than  the  airplane  frequency 

must  be  used  in  order  to  improve  the  damping.  The  possibility  of  using  • 
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a  negatively  (reverse)  geared  second-order  rate  autopilot  to  improve  the 
damping  seems  rather  surprising,  since  for  a  perfect  proportional  rate 
autopilot,  reverse  gearing  would  simply  decrease  the  effective  Cnr  and 

therefore  decrease  the  damping.  The  use  of  reverse  gearing  is  made 
possible  by  the  phase  relations  introduced  between  the  airplane  and 
rudder  motions  by  the  dynamics  of  the  second-order  autopilot.  Actually, 
this  possibility  is  no  more  surprising  than  the  fact  that  this  type  of 
autopilot  will  improve  the  damping  with  positive  gearing  only  when  the 
autopilot  natural  frequency  is  greater  than  the  airplane  frequency. 
Clearly,  this  restriction  arises  from  the  same  type  of  phase-relation 
requirement. 

The  use  of  a  negatively  geared  yaw  damper  would  seem  advantageous 
because  of  its  properties  in  a  steady  turn.  For  constant  yawing  velocity 
the  negatively  geared  yaw  damper  deflects  the  rudder  in  a  direction  to 
maintain  the  turn,  whereas  the  positively  geared  yaw  damper  must  be 
overridden,  either  by  the  pilot  or  by  the  boost  system.  However,  certain 
objections  to  the  use  of  negative  gearing  in  an  automatic  damper  for  use 
with  an  airplane  actually  make  such  use  impractical. 

The  main  objection  to  the  use  of  negative  gearing  is  based  on  the 
fact  that  the  equivalent  oscillator  represents  the  airplane  only  in  a 
given  flight  condition.  At  different  flight  conditions  the  character¬ 
istic  airplane  oscillation  has  different  values  of  damping  and  frequency; 
therefore,  the  airplane  is  represented  by  a  different  equivalent  oscil¬ 
lator  at  each  flight  condition.  To  design  an  automatic  damper  for  one 
flight  condition  only  is  impractical,  since  this  automatic  damper  may 
have  a  harmful  effect  on  the  damping  at  some  other  flight  condition. 
Autopilot  characteristics  must  therefore  be  obtained  by  some  compromise 
method  which  will  improve  any  practical  flight  condition.  Now,  it  can 
be  shown  that  the  regions  of  improved  damping  in  the  £cD0-plane  for  nega¬ 
tive  values  of  K  are  loops  resembling  a  reflection  in  the  A  =  0  axis 
of  the  unstable  loop  shown  in  figure  1.  These  loops  must  lie  in  a  rela¬ 
tively  narrow  range  of  cjd0  values,  since  they  are  confined  to  values 
of  0)0  lower  than  the  airplane  frequency.  Moreover,  for  a  given  magni¬ 
tude  of  gearing  ratio,  the  loop  for  any  damping  is  much  smaller  for  nega¬ 
tive  K  than  for  positive  K.  Figure  11  shows  that  the  loops  for  nega¬ 
tive  values  of  K  break  down  at  much  smaller  values  of  damping  than  the 
loops  for  positive  values  of  K  of  the  same  magnitude.  If  these  small 
loops  in  the  £o)0-plane  are  drawn  for  a  desired  amount  of  damping  for  two 
extreme  flight  conditions  with  different  natural  frequencies,  the  possi¬ 
bility  of  their  intersecting  in  a  region  of  the  £o>o-plane  which  would 
give  the  desired  damping  to  both  flight  conditions  is  relatively  small. 

For  positive  gearing,  on  the  other  hand,  improved  damping  can  be 
obtained  up  to  the  critical  damping  for  an  infinite  range  of  o>0  starting 
at  a  value  of'  <x>0  somewhat  greater  than  the  natural  airplane  frequency 
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in  the  particular  flight  condition.  There  is,  therefore,  an  infinite 
range  of  u>0  values  which  will  insure  at  least  the  critical  damping 
for  any  number  of  flight  conditions.  The  minimum  value  of  cjdd  neces¬ 
sary  is  somewhat  greater  than  the  highest  natural  frequency  of  any  of 
the  possible  flight  conditions  of  the  airplane. 


Maximum  damping  for  any  gearing  ratio.-  The  highest  maximum  damping 
and  the  corresponding  value  of  K  can  be  obtained  by  determining  where 
the  modes  corresponding  to  the  optimum  points  in  figures  10  and  11  become 

nonoscillatory.  This  condition  will  occur  for  P2  =  4Q.  By  using  this 
condition  in  equation  (17),  where  the  positive  sign  corresponds  to  posi¬ 
tive  K  and  the  negative  sign  to  negative  K,  the  following  expressions 
are  obtained: 


p- = 2K\F1 + x) 

p"“ = 2^°{\ Ffe  • x) 


(K  >  0) 


(K  <  0) 


(19) 


(20) 


Since  PQ  is  small  for  lightly  damped  airplanes,  the  positive  gearing 
ratio  gives  a  higher  value.  For  example,  for  PQ  =  0,  equation  (19) 
gives  a  value  approximately  six  times  as  large  as  equation  (20). 

Actually,  these  limiting  values  for  the  damping  are  of  only  academic 
interest  as  far  as  application  to  the  airplane  is  concerned,  since  they 
are  so  large  as  to  be  far  above  any  required  damping. 

Constant-damping  curves  in  Ko)0 -plane  with  damping  ratio  fixed.-  In 

order  to  obtain  a  more  complete  picture  of  the  effect  of  varying  gearing 
ratio  on  the  stability  of  the  system,  constant- damping  curves  were 
obtained  in  the  Kx>0-plane  with  £  fixed.  As  shown  in  appendix  A,  it 
is  necessary  to  solve  a  quadratic  equation  for  the  o>0  values,  which 
may  then  be  substituted  into  an  expression  for  K,  as  follows: 

(qo  -  R2  -  <jd2)o>02  -  2£u>0(P0  +  2R)(cd2  +  R2)  +  cn^  - 
(Qq  +  2PqR  +  2R2)(cd2  +  R2)  =  0  (21a) 
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K  =  +  po +  - 

Cl%  I 


(PQ  +  2R)a>02  +  2^0(q0  +  P0R  +  3R2) 


For  the  first  calculations,  5  =  0.3  was  assumed  as  a  reasonable 
value  of  damping  ratio  for  the  autopilot.  A  set  of  curves  was  obtained 
for  this  value  of  £ .  The  effect  of  varying  £  on  the  curves  in  the 
Kh)0-plane  was  then  investigated  by  obtaining  several  typical  curves  for 

5  =0.6  and  ^=0. 9* 

The  zero-damping  curve  for  5  =0.3  (shown  in  fig.  12)  indicates 
that  better-than-neutral  damping  can  be  obtained  for  any  value  of  cdq 
and  for  positive  or  negative  gearing.  This  is  true  for  any  damping  up 
to  the  airplane  damping.,  The  boundaries  for  less  than  the  airplane 
damping  are  of  no  practical  interest,  however,  and  figure  12  is  pre¬ 
sented  only  for  completeness. 

Figure  13  shows  the  airplane  damping  boundary,  which  is  a  simple, 
continuous  curve.  This  curve  alone  does  not  indicate  clearly  the  region 
that  defines  points  which  give  better  damping  than  that  for  the  air¬ 
plane  without  yaw  damper.  However,  the  axis  K  =  0  must  be  part  of  the 
boundary  also,  since  K  =  0  implies  no  autopilot,  which  means  that  the 
airplane  has  its  original  damping.  In  order  to  verify  that  the  region 
defined  by  the  hatching  in  figure  13  is  the  good  region,  a  curve  Was 
drawn  for  a  slightly  greater  damping  (^l/2  =  2.50  seconds).  This  curve 
is  shown  in  figure  l4  and  confirms  the  fact  that  the  region  insuring 
damping  greater  than  that  of  the  airplane  is  as  shown  in  figure  13 • 

The  results  of  figure  13  confirm  several  of  the  previous  conclu¬ 
sions  concerning  negative  gearing  which  were  obtained  from  figures  10 
and  11.  A  relatively  small  region  is  present  in  the  Kh>0-plane  in  which 
improved  damping  can  be  obtained  with  negative  K,  and  this  region  is 
confined  to  frequencies  less  than  the  airplane  frequency. 

Figure  l4  shows  that  the  regions  of  negative  K  which  will  give 
improved  damping  are  loops  in  the  Kh>0-plane.  Thus,  the  maximum  damping 
for  negative  K  may  be  obtained  from  the  cusp  point  corresponding  to 
the  breakdown  of  these  loops  for  any  values  of  5*  Since  5  is  con¬ 
stant,  a  particularly  simple  expression  can  be  obtained  for  Rmax.  ' 

If  5  and  ^o  are  use(i  instead  of  A  and  B,  equations  (l4)  and  (15) 
may  be  used  in  equation  (17)>  and  the  result  is 


+  k(2Qq  +  3PdR  +.  ^R2) 


(21b) 
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%  = 


1  -  £ 


(K  >  0)  (22) 


^5  +  ^ 
1  +  £ 


(K  <  0)  (23) 


for  the  values  at  the  cusp  points.  These  values  are  given  by  the 
two  curves  in  figure  10.  The  maximum  damping  for  negative  K  'at  any 
value  of  £  may  then  be  obtained  by  using  equation  (23)  in  equa¬ 
tion  (l4)  to  yield 


(K  <  0) 


If  the  values  of  PQ  and  Qq  for  the  Dutch  roll  of  the  airplane  under 
consideration  are  used  and  £  is  taken  equal  to  0.3,  equation  (24) 
shows  that  the  maximum  damping  with  negative  gearing  corresponds  to 
T1/2  =  1.0  second. 

Figure  15  shows  the  curve  for  this  value  of  damping.  The  cusp 
point  for  negative  gearing  occurs  at  0>o  =  3.96  radians  per  second  and 

K  =  -0.035  deg/deg/sec.  For  positive  gearing,  on  the  other  hand,  an 
infinite  range  of  o)0  values  which  will  give  better  than  this  damping 
with  £  =  0.3  is  seen  to  be  available.  This  bears  out  the  previous 
statement  that  the  design,  with  regard  to  damping  ratio  and  natural  fre¬ 
quency,  of  a  compromise  autopilot  which  will  improve  a  variety  of  flight 
conditions  is  less  restricted  when  positive  gearing  is  used. 

Since  the  curve  of  figure  15  is  typical  of  the  curves  for  dampings 
somewhat  greater  than  the  airplane  damping,  it  will  be  discussed  in  more 
detail.  The  curve  has  two  separate  branches  and  the  discontinuity  occurs 
between  m  =  4.8  radians  per  second  and  cd  =  5  radians  per  second.  This 

fr'i+.'iml  ■FT*pnii»nr"u-  op+nollir  io  —  v/n  _  w2  /  ~~~  /  oi  ^  \  m 


critical  frequency  actually  is  (Ocr  -  \/Q0  -  R^  (see  eq.  (21a)).  From 
what  has  been  said  previously,  the  portion  of  the  curve  of  greatest 
interest  is  the  second  half  (for  positive  K) .  This  starts  at  the  criti- 

J>  +  2R 

cal  frequency  at  u>0  =  °°  and  K  =  -— - .  This  value  of  K  is  the 

Cl 

value  for  which  an  ideal  autopilot  (one  with  no  lags)  would  yield  the 
given  damping,  and  will  be  denoted  as  K0.  The  bottom  part  of  the  bound¬ 


ary  can  be  seen  to  represent  the  values  of 


and  cnn  which  make  the 
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airplane  mode  have  the  given  damping,  since  the  critical  frequency  is 
near  the  airplane  frequency.  On  the  other  hand,  the  top  part  of  the 
boundary  represents  the  values  of  K  and  o)0  at  which  the  autopilot 
mode  has  the  given  damping,  since  the  values  of  u>  along  this  part  of 
the  curve  are  very  close  to  the  corresponding  values  of  <*>0. 

Figure  l6  shows  the  boundary  for  Tj_/2  =0*^0  second.  The  loops 
for  negative  damping  have  disappeared.  Thi,s  curve  is  the  cusp  curve 
for  positive  K,  obtained  by  using  equation  (22)  in  equation  (14) .  This 
fact  has  little  practical  significance,  however,  since  there  is  no  funda¬ 
mental  change  in  the  shape  of  the  curves  at  this  damping. 

For  large  values  of  damping,  further  changes  occur  in  the  type  of 
damping  curves  in  the  Rn0- plane .  For  small  values  of  R,  equation  (21a) 
gives  real  solutions  for  cnQ  with  any  value  of  cn.  However,  for  larger 
values  of  R  the  discriminant  of  equation  (21a)  changes  sign  for 'certain 
combinations  of  o>  and  R,  so  that  no  real  solutions  for  0>o  exist. 
Setting  this  discriminant  equal  to  zero  yields 

U2  +  R2)2  +  [?(Po  +  2R)2  "  2R(Po  +  2R)  -  2Q^]  U2  +  R2)  + 

QojQc,  +  2R(P0  +  2R)]  =  0  (25) 

2  2 

When  this  equation  has  a  positive,  real  root  for  the  quantity  u>  +  R  , 
this  root  determines  the  range  of  real  values  of  u>  for  which  the  solu¬ 
tions  for  (So'  in  equation  (21a)  are  complex.  The  value  of  R  for 
which  real  roots  occur  in  equation  (25).  is  obtained  by  equating  the 
discriminant  of  this  equation  to  zero  to  give 

(1  -  ?V  -  (1  -  ♦  ?2(|  po¥  -  «o)  -  .0  <26> 

For  the  airplane  being  considered,  with  =  0.3,  the  values  of  R 
obtained  from  equation  (26)  are  R  =  -1.58  and  R  =  1.6l.  Since  only 
the  positive-damping  boundaries  are  of  interest,  only  the  negative  value 
(which  corresponds  to  Tp/2  =0.44  second^  need  be  considered.  For  this 

value  of  damping,  the  range  of  values  given  by  equation  (25)  becomes  a 
single  value.  For  greater  damping,  equation  (25)  gives  a  finite  range 
of  values  of  o>  which  will  not  occur  at  any  real  value  of  o>0.  Thus, 

no  value  of  co  in  this  range  can  occur  for  the  required  damping  when 
£  =  0.3.  The  constant-damping  boundaries  for  these  larger  values  of 
damping,  which  are  characterized  by  the  absence  of  a  given  range  of 
a)  values,  are  of  a  somewhat  different  type  from  the  curves  for  lower 

damping. 
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As  an  example  of  the  type  of  curve  which  occurs  for  large  damping 
values,  the  boundary  for  T^/2  =  0*25  second  (R  =  -2.77)  is  presented 

in  figure  17-  Substituting  £  =  0.3  and  R  =  -2.77  into  equation  (25) 
gives  co  =  4.3  radians  per  second  and  co  =  6.3  radians  per  second. 

That  is,  no  values  of  co  in  the  range  4.3  <  co  <  6.3  radians  per  second 
occur  on  this  curve.  Some  of  the  significant  features  of  this  curve 
will  now  be'  described,  since  it  is  typical  of  the  higher- damping  curves. 

The  significant  portion  of  this  curve  defines  the  same  general  type 
of  wedge-shaped  good  region  as  was  present  for  lower  dampings.  On  the 
lower  part  of  the  significant  boundary,  the  frequency  of  the  mode  with 
the  given  damping  is  in  the  range  4.0  <  cjd  <  4.3  radians  per  second; 
that  is,  the  frequency  varies  between  the  critical  frequency  and  the 
smaller  frequency  given  by  equation  (25).  These  frequencies  clearly 
represent  the  airplane  mode.  At  the  point  of  the  wedge  the  previously 
discussed  discontinuity  appears.  At  this  point  both  the  airplane  and 
autopilot  modes  have  the  same  damping.  On  the  upper  part  of  the  bound¬ 
ary  the  autopilot  mode  has  the  required  damping  and  the  airplane  mode 
has  higher  damping. 

Figure  18  is  a  collection  of  the  significant  portions  of  the  damping 
curves  previously  discussed.  This  figure  shows  that,  as  the  required 
damping  increases,  the  wedge-shaped  region  in  the  positive  Km0-plane 
which  insures  this  amount  of  damping  moves  upward  and  to  the  right.  The 
figure  also  shows  that  the  damping  obtainable  with  a  second-order  yaw 
autopilot  cannot  always  be  increased  merely  by  increasing  the  gearing 
ratio.  The  reason  is  that,  fqr  a  given  value  of  cn0,  increasing  K 

beyond  a  certain  value  makes  the  autopilot  mode  less  stable.  Clearly, 
more  damping  can  be  obtained  by  increasing  o>0  at  the  same  time  that 
K  is  increased.  These  results  confirm  the  statement  made  in  the  dis¬ 
cussion  of  figure  10  that  it  may  be  necessary  to  increase  the  natural 
frequency  when  the  gain  is  increased. 

It  is  important  to  remember  that,  as  in  the  £<n0-plane,  two  oscil¬ 
lations  are  present  in  the  regions  of  most  interest  in  the  Kh>0-plane. 

The  typical  frequencies  given  on  the  individual  boundaries  indicate  that 
the  portions  of  the  wedge-shaped  boundaries  where  K  is  high  correspond 
to  the  autopilot  mode  having  the  given  damping,  whereas  the  parts  where 
K  is  low  correspond  to  the  airplane  mode  having  the  given  damping.  . 

Figure  18  seems  to  indicate  that  an  infinite  amount  of  damping 
might  be  obtainable  by  simultaneously  increasing  K  and  cjdg,  in  contra¬ 
diction  to  the  discussion  concerning  figures  10  and  11.  However,  further 
changes  occur  in  the  type  of  curve  at  larger  dampings.  Equation  (21a) 
indicates  that  a  change  might  be  expected  when  R  >  since  the  quan- 

2  2 

tity  Q0  -  R  -  u>  does  not  go  through  zero  for  any  value  of  cd.  That 
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is,  there  no  longer  is  a  critical  frequency  o>cr  =  V  Qo  -  R2-  Moreover, 
the  good  regions  as  previously  shown  do  not  take  into  account  the  bound¬ 
ary  of  equal  roots.  That  is,  within  the  good  regions  no  oscillation  has 
less  damping  than  the  indicated  amount;  however,  a  real  root  with  less 
damping  than  the  indicated  amount  may  he  present  at  points  above  the 
boundary  of  equal  roots.  The  explanation  is  as  follows.  At  the  <x>  =  0  ■ 
point  of  a  given  high-damping  curve  (which  point  is  on- the  boundary  of 
equal  roots)  there  are  two  equal  real  rpots  with  the  indicated  amount 
of  damping.  However,  at  points  above  the  boundary  of  equal  roots,  one 
of  these  real  modes  will  decrease  in  damping.  The  oscillatory  mode 
which  breaks  down  into  two  real  modes  at  the  boundary  of  equal  roots  in 
the  Kb0-plane  is  the  airplane  mode,  since,  as  can  be  seen  from  fig¬ 
ure  18,  oscillations  at  autopilot  frequency  do  occur  above  the  boundary 
of  equal  roots.  Thus,  for  points  above  this  boundary  the  airplane  mode 
becomes  nonoscillatory.  -The  physical  reason  for  this  phenomenon  is  that 
large  gearing  ratios  cause  the  airplane  mode  to  be  overdamped.  The 
boundary  of  equal  roots  in  the  £o3Q-plane,  on  the  other  hand,  corresponds 
to  large  £  values.  Therefore,  as  can  be  seen  from  the  curves,  the 
mode  which  is  overdamped  at  points  above  this  boundary  is  the  autopilot 
mode.  (  > 

Strictly  speaking,  the  boundary  of  equal  roots  should  be  considered 
as  the  upper  limit  of  the  good  regions  shown  in  figure  18.  Actually, 
this  restriction  is  necessary  only  for  the  regions  of  very  high  damping. 

For  example,  when  the  point  K  =  0.60  deg/ deg/sec  and  03o  =  21.5  radians 

per  second  is  taken  on  the  boundary  for  T ]_/£  =  0.6o  second  in  fig¬ 
ure  18  and  the  characteristic  roots  are  found,  the  autopilot  mode  has 
the  given  damping  T^yp  =  0*^0  second  at  o>  =  21  radians  per  second, 

whereas  the  airplane  oscillation  breaks  down  into  two  well-damped  non¬ 
oscillatory  modes  with  T -^/q  ~  0.22  second  and  T \j2.  =  0.09  second. 

Thus,  for  this  moderate  damping,  points  -above  the  boundary  of  equal 
roots  still  give  the  required  damping. 

For  extremely  large  values  of  damping,  however,  the  boundary  of 
equal  roots  becomes  important.  In  fact,  the  maximum  damping  of  the 
complete  system  at  a  fixed  value  of  t,  is  the  largest  negative  value 
of  R  occurring  on  the  boundary  of  equal  roots.  The  points  on  the 
boundary  of  equal  roots  are  obtained  by  using  cn  =  0  in -equations  (21a) 
and  (21b).  The  largest  value  of  R  occurring  on  the  boundary  is 
obtained  by  setting  the  discriminant  of  equation  (21a)  equal  to  zero 
with  a)  =  0.  This  procedure  gives  the  quartic  in  R  (from  eq.  (25) 
with  (jd  =  0) : 

{hf  -  3)r4  +  2P0(2£2  -  l)R3  +  (?2P02  +  2Qo)r2  +  2P0QoR  +  Q02  =  0  (27) 
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The  largest  negative  real  root  of  this  equation  is  the  value  of  the 
maximum  damping.  For  the  example  considered,  this  value  is  R  =  -5-11 
and  the  corresponding  Tp/p  =  0-1^  second.  For  this  value  of  damping, 

the  good  region  becomes  a  single  point  on  the  boundary  of  equal  roots. 

The  airplane  mode  is  a  double  real  root  equal  to  -5-11,  and  the  auto¬ 
pilot  mode  is  an  oscillation  with  the  same  amount  of  damping.  The 
autopilot  'characteristics  for  this  maximum-damping  point  were  used  in 
calculating  the  motion  subsequent  to  a  5°  sideslip  disturbance  (fig.  19). 
For  practical  purposes,  the  airplane  motion  can  be  considered  as  non- 
oscillatory,  since  the  autopilot  oscillator  is  of  too  high  a  frequency 
to  have  an  appreciable  effect  on  the  airplane.  The  well-damped  nature 
of  the  Dutch  roll  motion  (which  has  now  become  nonoscillatory)  is  evi¬ 
dent  in  the  (3-motion.  The  slow  return  of  the  roll  motion  is  due  to  the 
spiral  mode.  Although  the  lightly  damped  spiral  mode  is  generally  not 
considered  troublesome,  it  is  necessary  to  keep  in  mind  that  the  dis¬ 
cussion  in  this  paper  deals  only  with  the  improvement  of  the  damping  of 
the  Dutch  roll  mode,  and  that  the  two  aperiodic  modes  in  the  airplane's 
lateral  motion  have  been  ignored. 

Effect  of  varying  £  on  curves  in  the  Kh)0 -plane.-  In  order  to 
obtain  an  idea  of  the  effect  of  changing  £  on  the  curves  in  .the 
Ko)0-plane,  a  comparison  of  the  boundary  of  equal  roots,  the  airplane 
damping  boundary,  and  the  boundary  for  T =  0.6o  second  is  presented 

in  figure  20  for  £  =  0.3,  0.6,  and  0.9.  Although  the  region  of  improved 
damping  with  negative  gearing  increases  in  size,  the  narrowness  of  the 
frequency  range  and  the  other  difficulties  previously  mentioned  still 
make  the  use  of  negative  gearing  impractical.  The  increased  slope  of 
the  upper  part  of  the  boundaries  simply  implies,  as  would  be  expected, 
that  for  larger  values  of  £,  larger  values  of  K  are  required  At  any 
value  of  (J0o  to  make  the  autopilot  mode  become  unstable.  The  variation 
in  the  position  of  the  boundary  of  equal  roots  simply’  implies  that  the 
more  highly  damped  mode  associated  with  each  point  on  a  given  boundary 
tends  to  become  critically  damped  at  lower  gearings  as  £  increases. 

None  of  these  variations  appear  to  be  very  important  practically. 
Indeed,  the  most  important  fact  about  the  effect  of  varying  £  is  that 
the  bottom  part  of  the  boundary  approaches  the  same  value  of  K  at 
large  values  of  cn0  for  all  values  of  £.  This  behavior  is  due  to  the 
fact  that  the  asymptotic  value  of  K  is  K0,  the  value  required  to 
obtain  the  given  amount  of  damping  with  an  ideal  autopilot.  Now,  the 
bottom  part  of  the  boundary  can  be  seen  to  be  the  important  part,  since 
it  gives  the  lowest  value  of  K  for  which  the  system  has  the  required 
damping;  also,  the  mode  which  attains  the  required  damping  at  'the  points 
on  this  part  of  the  boundary  is  the  airplane  mode.  The  flatness  and 
invariance  with  £  of  the  bottom  part  of  the  boundary  imply  that  the 
minimum  value  of  K  for  which  a  given  damping  may  be  attained  is  rela¬ 
tively  invariant  for  changes  in  t,  and  o>0  at  fairly  large  values 
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of  o)0.  This  fact  will  he  used  in  the  next  section  to  obtain  a  simple 
rule-of-thumb  method  for  designing  an  efficient  yaw-rate  damper  which 
will  insure  a  required  damping  for  several  extreme  flight  conditions  of 
an  airplane  using  low  gearing. 

Design  of  Compromise  Yaw-Rate  Damper 


.  For  the  sake  of  simplicity,  the  autopilot  to  be  designed  will  be 
assumed  to  provide  only  a  yawing  moment  proportional  to  yawing  velocity, 
as  in  the  previous  analysis.  It  may  therefore  be  called  a  Cnr  auto¬ 
pilot.  Values  of  K,  £,  and  oo0  which  will  efficiently  improve  the 
damping  of  an  airplane  in  various  flight  conditions  are  desired.  The 
criterion  of  the  efficiency  of  the  autopilot  will  be  that  the  gearing 
ratio  required  for  the  autopilot  must  be  small.  This  means  that  the 
autopilot  power  required  will  be  small.  Moreover,  the  use  of  small 
gearing  makes  the  yaw  damper  easier  to  override  in  steady  turns. 

Table  II  gives  the  parameters  used  for  the  three  flight  conditions 
of  the  airplane  chosen  for  the  present  example.  Case  I  is  a  high-lift- 
coefficient,  low -wing-loading  landing  condition  at  sea  level.  Case  II 
is  a  low-lift-coefficient,  medium-wing-loading  cruising  condition  at 
30,000  feet.  To  complete  the  picture,  case  III  is  a  high-lift- 
coefficient,  high-wing- loading  cruising  condition  at  30,000  feet.  These 
cases  will  serve  as  examples  to  illustrate  the  method. 

Table  II  shows  that  the  Dutch  roll  oscillation  in  case  II  is  very 
poorly  damped,  since  it  requires  almost  7  seconds  to  damp  to  one-half 
amplitude.  Although  the  other  two  cases  are  not  so  poorly  damped,  they 
are  still  unsatisfactory.  Calculated  motions  for  the  three  cases  in 
response  to  a  5°  disturbance  in  sideslip  are  shown  in  figure  21.  No 
attempt  will  be  made  to  set  up  any  complicated  criteria  for  adequate 
damping.  Instead,  the  criterion  chosen  for  purposes  of  illustration 
will.be  that  the  Dutch  roll  oscillation  should  damp  to  one-half  ampli¬ 
tude  in  1  second  or  less  at  any  flight  condition.  Actually,  the  auto¬ 
pilot  may  be  designed  to  insure  a  different  amount  of  damping  for  each 
flight  condition,  in  case  one  of  the  flight  conditions  is  required  to 
be  more  stable  than  another. 


The  fundamental  problem  is  to  find  the  set  of  points  (K,<dq,  £)  which 
satisfy  the  given  damping  criterion  for  all  three  flight  conditions. 
Moreover,  the  minimum  possible  gearing  is  desired.  The  value  of  K0 
for  T^/2  =  1  second  is  calculated  for  each  flight  condition.  From 
the  equivalent-oscillator  analysis,  the  value  of  K0  is  shown  to  be 
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Substitution  from  table  II  gives  the  values  of  K0  as  0.1213,  O.0698, 

/■ 

and  0.137^  deg/deg/sec  for  cases  I,  II,  and  III,  respectively. 

If  the  largest  of  these  values  of  K0  (corresponding  to  case  III) 
is  chosen  for  the  gearing  of  the  autopilot,  the  curve  for  Tp/p  =  1  second 
in  the  £<jo0-plane  will  be  the  critical  curve  for  case  III,  but  for  cases  I 
and  II  the  curve  for  T p/p  =  1  second  will  be  a  curve  between  the  air¬ 
plane  damping  and  the  critical  damping.  Each  of  these  three  curves 
defines  an  infinite  region  in  the-  upper  right  portion  of  the  £o>0-plane 
as  the  region  insuring  better  damping.  Therefore,  in  the  £u>0-plane  an 
infinite  region  of  points  which  are  common  to  these  three  regions  must' 
exist  and  can  be  used  with  this  value  of  K  to  obtain  better  damping 
than  T-^2  =  1  second  for  all  three  flight  conditions.  A  plot  of  these 

regions  is  shown  in  figure  22(a),  where  K  =  0.l4  deg/deg/sec  has  been 
used.  The  plot  shows  that  any  set  of  values  in  the  £(i>0 -plane  which 
insures  Tj/2  =  ^  second  for  case  II,  the  high-frequency  case,  will 

also  insure  this  amount  of  damping  for  the  other  two  cases.  However, 
any  point  in  this  region  is  at  a  value  of  cDq  considerably  above  the 
optimum  for  cases  I  and  III.  As  was  previously  pointed  out,  for  such 
values  of  oq0  the  minimum  value  of  K  for  which  a  given  damping  may 
be  obtained  does  not  differ  much  from  K0.  The  implication  is  that  the 
minimum  value  of  K  for  which  Tp/2  =  1  second  may  be  insured  for 
case  III,  while  at  the  same  time  Tp/2  remains  less  than  1  second  for 

case  II,  is  not  much-  less  than  0.137^-  deg/ deg/sec.  To  confirm  this 
hypothesis,  the  necessary  points  on  the  curves  in  the  £o>0-plane  for 
Tf/2  =  1  second  were  obtained  for  K  =  0.12  deg/deg/sec  and 

K  =  0.10  deg/deg/sec,  and  the  significant  parts  of  these  curves  are 
shown  in  figures  22(b)  and  22(c). 

For  values  of  K  <0.137^  deg/deg/sec  the  curve  for  Tpy^  =  1  second 
for  case  III  must  be  a  loop,  since  the  damping  is  higher  than  the  critical 
damping.  Figure  22(b)  shows  that  when  K  =  0.12  deg/deg/sec  this  loop 
is  still  large  enough  to  intersect  the  region  where  Tp/g  <  1  second  for 

case  II.  The  hatched  region  on  figure  22(b)  is  the  region  which  will 
insure  Tp <  1  second  for  all  three  flight  conditions  at  K  =  0.12 

deg/deg/sec.  Figure  22(c)  shows  that  the  loop  for  case  III  has  become 
too  small' to  intersect  the  region  where  Tp <  1  second  for  case  II 

when  K  =  0.10  deg/deg/sec.  Therefore,  a  value  of  Tp <  1  second 
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cannot  be  obtained  for  cases  II  and  III  simultaneously  when  K  =  0.10 
deg/deg/sec.  The  curve  for  case  I  is  not  shown  in  this  figure  since , 
as  seen  in  figures  22(a)  and  22(b) ,  it  is  not  necessary. 


The  discussion  of  figure  22  shows  that,  in  order  to  improve  all 
the  flight  conditions  in  the  present  example  with  a  single  second-order 
autopilot,  a  gearing  ratio  has  to  be  used  which  is  almost  as  large  as 
would  be  needed  for  an  ideal  autopilot  that  would  stabilize  all'  the 
flight  conditions.  This  difficulty  arises  because  the  low-frequency 
conditions  require  the  highest  values  of  K0,  as  can  be  seen  by 
examining  the  approximate  expression  for  the  frequency  of  the  Dutch 
roll  oscillation  at  any  flight  condition  and  the  associated  expression 
for  K0 : 


Because  of  the  possibility  of  variations  in  C^,  Cng,  an<i  P  "  Po 

in  the  various  flight  conditions,  it  is  not  possible  to  state  that  in 
general  the  low-frequency  flight  conditions  will  require  the  higher 
values  of  K0,  as  is  true  in  the  present  example.  When  the  highest 
required  value  of  K0  occurs  at  the  high-frequency  conditions,  the  ' 
minimum  value  of  K  for  the  compromise  autopilot  may  be  considerably 


below  this  highest  required  value,  because  the  region  of  overlap  similar 
to  that  shown  in  figure  22(b)  will  include  the  optimum  point  for  the 


high-frequency  condition. 


In  any  case,  the  characteristics  of  the  compromise  autopilot  can  . 
be  obtained  by  plotting  the  regions  insuring  the  required  damping  for 
the  various  conditions  at  various  values  of  K  below  the  maximum  value 
of  K0,  as  is  done  in  figure  22.  The .value  of  K  which  leads  to  a 
small  overlap  region  is  then  the  minimum  compromise  gearing,  and  the 
values  in  the  overlap  region  of  the  £to0-plane  define  the  possible  values 
'for  £  and  co0  of  the  compromise  autopilot. 
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A  convenient  method  is  available  for  determining  the  range  of 
values  of  K  to  be  used  in  obtaining  the  overlap  region.  The  values  of 
K0  and  the  optimum  value  of  K  which  will  yield  the  required  damping 
for  each  of  the  flight  conditions  are  calculated.  Then  the  compromise 
value  of  K  must  lie  between  the  highest  optimum  value  of  K  and  the 
highest  value  of  K0.  In  making  the  plots  it  may  soon  become  clear  that 
some  of  the  flight  conditions  do  not  need  to  be  considered,  as  happened 
with  case  I  in  the  present  example. 


In  the  present  example,  it  has  been  shown  that  a  gain  much  smaller 
than  the  largest  value  of  K0  necessary  cannot  be  used,  unless  a  vari¬ 
able  gain  is  available.  If  a  constant  gain  must  be  used  at  all  flight 
conditions  of  such  an  airplane,  a  simple  rule  of  thumb  for  designing 
the  autopilot  would  be  to  choose  the  largest  value  of  K0  (corre¬ 
sponding  to  the  low-frequency  conditions)  and  values  of  £  and  cn0 
near  the  optimum  point  of  the  high-frequency  condition.  In  this  way 
some  advantage  is  derived  from  the  second-order  characteristics  of  the 
autopilot  in  that  the  damping  obtained  from  the  autopilot  is  much  larger 
for  the  high-frequency  conditions  than  that  which  would  result  from  an 
ideal  autopilot  with  this  gearing. 


In  applying  this  simple  method,  0.l4  deg/deg/sec  was  chosen  for 
the  value  of  K,  and  the  optimum  point  was  found  for  case  II.  The  opti¬ 
mum  point  was  t,  =  0.523  and  o>0  =  9*^9  radians  per  second.  Figure  23 
shows  the  calculated  motions  for  the  three  flight  conditions  with  this 
autopilot,  subsequent  to  a  5°  sideslip  disturbance.  A  comparison  of 
figure  23  with  figure  21  reveals  that  the  stability  of  all  three  flight 
conditions  is  improved.  In  particular,  case  II  is  greatly  improved, 
while  cases  I  and  III  both  have  T^/g  slightly  less  than  1  second. 


Since  the  results  of  figure  2l(b)  show  that  a  gearing  of 
0.12  deg^deg/sec  could  have  been  used,  the  optimum  point  for  case  II 

was  calculated  for  this  gearing  and  found  to  be  £  =  0.485  and 
cjo0  =  8.8l  radians  per  second.  Figure  2b  shows  the  motions  in  the  three 
flight  conditions  with  these  autopilot  characteristics.  A  comparison 
of  this  figure  with  figure  23  reveals  that  the  damping  obtained  with 
this  autopilot  is  only  slightly  less  than  the  damping  obtained  when 
K  =  0.14  deg/deg/  sec. 


If  the  low-frequency  conditions  require  little  improvement,  the 
problem  approaches  that  of  improving  only  a  single  flight  condition, 
so  that  the  optimum-point  characteristics  may  be  used  to  decrease  the 
gearing  necessary.  For  example,  suppose  that  T-^/g  =  1.5  seconds  had 


been  considered  satisfactory  for  the  low-frequency  conditions.  The 
value  of  K0  necessary  to  obtain  this  damping  is  0.075  degy/deg/sec, 
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and  the  corresponding  optimum  point  for  case  II  is  £  =0.389  and 
<jo0  =  7> ^*0  radians  per  second.  The  improvement  of  the  stability  of  the 
various  flight  conditions  may  be  seen  from  the  motions  presented  in  fig- 
tire  25-  In  these  motions  the  value  of  T^/2  ^or  cases  I  and  III  is 

between  1  and  1.5  seconds,  but  the  value  of  Tj/g  for  case  II  is  well 
under  1  second. 

It  is  clear  from  the  previous  discussion  that  for  the  present 
example  an  autopilot  with  fixed  characteristics  would  require  a  large 
gain  (corresponding  to  the  largest  value  of  K0)  in  order  to  improve 
both  landing  and  cruising  flight  conditions.  4A  variable  gain  would 
therefore  be  desirable,  so  that  this  excessive  gearing  would  not  have 
to  be  used  at  the  higher  speeds .  Although  it  would  be  impractical  to 
expect  that  all  three  autopilot  parameters  should  be  variable  in  flight, 
as  would  be  desirable  to  obtain  the  optimum  autopilot  for  each  flight 
condition,  it  would  be  relatively  simple  to  make  the  gain  variable .  The 
value  of  £  chosen  would  be  that  of  the  optimum  point  for  the  high- 
frequency  condition  (as  calculated  for  the  low-gearing  value),  and  the 
value  of  u>0  would  be  preferably  slightly  above  the  optimum  value  of 

o>0  for  the  high-frequency  condition.  (Because  of  the  rapid  change  in 
damping  at  values  of  cnQ  less  than  the  optimum,  in  practice  a  safety 
factor  should  generally  be  added  to  this  value.)  In  such  cases,  prob¬ 
ably  only  two  gain  positions  would  be  necessary  -  a  high  gain  for  the 
low-velocity  conditions  and  a  low  gain  for  the  high-velocity  conditions. 
Thus,  low  gain  could  be  used  in  the  high-velocity  flight  conditions,  so 
that  the  adverse  effects  of  the  autopilot  in  steady  turns  at  high 
velocity  would  be  small. 


Validity  of  Assuming  Pure  Cn^  Autopilot 

A  few  final  remarks  will  be  made  concerning  the  assumption  that  the 
autopilot  is  sensitive  to  yawing  velocities  only,  and  provides  yawing 
moments  only.  For  practical  autopilot  installations,  the  sensing  device 
(usually  a  gyroscope)  is  fixed  in  the  airplane.  The  device  is  therefore 
sensitive  to  yawing  velocities  about  some  axis  fixed  in  the  airplane. 

The  equations  of  motion,  however,  are  set  up  with  respect  to  the  sta¬ 
bility  axes.  For  various  flight  conditions,  the  angle  of  inclination 
between  the  gyro  axes  and  the  stability  axes  varies  because  of  the 
varying  angle  of  attack.  This  angle  of  inclination  makes  the  gyroscope 
sensitive  also  to  rolling  velocity  about  the  stability  axes.  In  addi¬ 
tion,  the  displacement  of  the  autopilot-actuated  control  surface  from 
the  longitudinal  stability  axis  gives  rise  to  rolling  moments.  As  shown 


in  reference  10,  these  two  effects  cause  increments  in  C 
C  2^  due  to  the  autopilot,  in  addition  to  the  expected 


np> 


"n  T 


C ir,  and 
increment . 


The  Cn  effect  is  the  most  important  of  these  in  affecting  the  Dutch 


roll  stability  of  the  airplane,  and  this  effect  will  now  be  discussed. 
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When  the  stability  of  a  single  flight  condition  is  to  he  improved, 
the  Cn  effect  may  he  removed  hy  orienting  the  gyro  axis  along  the 

ir  I  ■ 

flight  path,  and  no  modification  of  the  results  is  necessary.  If  an 
angle  of  gyro  inclination  does  exist,  the  general  method  of  setting  up 
the  constant-damping  boundaries  is,  of  course,  still  valid.  In  order 
to  take  account  of  the  Cn  effect,  Dty  is  simply  replaced  in  equa¬ 
tion  (6)  hy  the  proper  linear  combination  of  Dty  and  D0,  as  shown  in 
reference  10.  Because  of  the  presence  of  the  D term  all  three  degrees 
of  freedom  must  he  used,  and  the  equivalent-oscillator  simplification  can 
no  longer  he  used.  Thus,  the  method  is  more  complicated  when  the  Cn 

P 

effect  of  the  autopilot  is  not  small  enough  to  he  neglected.  When  the 
Cnp  effect  becomes  excessively  large,  the  problem  becomes  even  more 

complicated  because  of  the  fact  that  the  damping  of  the  aperiodic  modes 
becomes  important,  and  curves  of  constant  damping  for  these  modes  (corre¬ 
sponding  to  real  roots)  have  to  be  plotted.  In  fact,  the  aperiodic  modes 
may  combine  to  form  another  oscillation. 

In  considering  the  problem  of  simultaneously  stabilizing  various 
flight  conditions  which  have  various  angles  of  attack,  the  Cn  effect 

must  always  be  considered  unless  the  gyro  axis  can  somehow  be  rotated 
so  that  it  is  always  parallel  to  the  relative-wind  axis.  The  Cnp  effect 

can  sometimes  be  ignored,  if  the  angle -of- attack  range  of  the  flight  con¬ 
ditions  is  small,  by  choosing  an  orientiation  of  the  gyroscope  in  such  a 
direction  that  all  the  angles  of  inclination  are  small.  Because  of  the 
extreme  complexity  of  the  problem  when  the  Cnp  effects  must  be  con¬ 
sidered,  this  paper  is  confined  to  the  consideration  of  Cnr  effects 
only. 


CONCLUDING  REMARKS 

The  damping  of  an  oscillatory  system  that  makes  use  of  a  second- 
order  rate  damper  with  a  given  gearing  ratio  K  can  be  improved  by 
adjusting  the  shape  of  the  autopilot  frequency  response  (that  is,  the 
damping  ratio  and  the  natural  frequency  of  the  autopilot) .  For  the 
purpose  of  determining  the  effect  of  a  second-order  yaw-rate  damper  on 
the  damping  of  the  Dutch  roll  motion  of  an  airplane  in  a  given  flight 
condition,  the  airplane  may  be  represented  as  an  equivalent  oscillator. 

By  using  this  equivalent-oscillator  concept,  the  optimum  shape  of  the 
autopilot  frequency  response  corresponding  to  a  given  gearing  ratio  or 
required  damping  may  be  obtained  from  a  simple  set  of  equations  for  any 
flight  condition.  The  gearing  ratio  necessary  to  obtain  a  given  amount 
of  damping  when  the  damping  ratio  and  natural  frequency  of  the  autopilot 
are  hear  their  optimum  values  is  considerably  less  than  the  gearing  ratio 
necessary  to  obtain  the  same  amount  of  damping  with  an  idealized  (no-lag) 
autopilot . 
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The  problem  of  designing  a  second-order  yaw  damper  for  an  airplane 
which  requires  improvement  in  damping  for  several  flight  conditions  is 
more  complicated,  since  for  each  flight  condition  the  airplane  is  repre¬ 
sented  by  a  different  equivalent  oscillator.  However,  a  simple  method 
of  compromise  is  derived  for  flight  conditions  in  which  the  effect  of 
the  yaw-damper  sensitivity  to  rate  of  roll  is  small  enough  to  be 
neglected. 

The  effects  of  a  second-order  yaw  damper  on  the  stability  of  any 
given  flight  condition  of  an  airplane  can  be  obtained  by  examining  the 
constant-damping  curves  in  the  plane  of  the  dampifig  ratio  and  natural 
frequency  of  the  autopilot  and  in  the  plane  of  the  gearing  ratio  and 
natural  frequency  of  the  autopilot.  Theoretically,  any  given  flight 
condition  may  be  stabilized  by  using  either  positive  or  negative  gearing. 
When  negative  gearing  is  used,  the  autopilot  natural  frequency  must  be 
less  than  the  airplane  frequency.  The  use  of  negative  gearing  is  shown 
to  be  impractical,  however.  For  positive  gearing,  the  autopilot  natural 
frequency  must  be  greater  than  the  airplane  frequency . 

I 

1  For  fixed  positive  gearing,  there  is  an  infinite  number  of  combi¬ 
nations  of  autopilot  natural  frequency  and  damping  ratio  for  which  the 
second-order  autopilot  gives  better  damping  than  an  ideal  autopilot  of 
the  same  gearing.  For. fixed  positive  damping  ratio  of  the  autopilot, 
there  is  a  range  of  values  of  positive  gearing  ratio  which  will  provide- 
a  given  damping  to  the  system  at  any  autopilot  natural  frequency  greater 
than  the  airplane  frequency.  Increasing  the  gearing  ratio  of  the  auto¬ 
pilot  to  excessive  values  will  always  cause  the  autopilot  mode  of  oscil¬ 
lation  to  become  unstable  for  a  given  damping  ratio  and  natural  frequency 
of  the  autopilot.  If  larger  gearing  ratios  are  to  be  used  in  order  to 
obtain  higher  stability  for  the  airplane  mode,  the  damping  ratio  or  the 
natural  frequency  of  the  autopilot,  or  both,  must  be  increased.  Expres¬ 
sions  are  derived  for  the  maximum  damping  under  various  conditions. 

Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 

Langley  Field,  Va.,  October  2,  1952. 
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APPENDIX  A 


EXPRESSIONS  FOR  CONSTANT -DAMPING  CURVES  AND  DISCUSSION 
OF  EQUIVALENT -OSCILLATOR  CONCEPT 


Since  the  characteristic  equation  of  the  airplane-damper  system  is 
of  sixth  degree,  as  given  in  equation  (7),  it  may  he  written 

F(D)  =  D6  +  A5D5  +  A^D4  +  A3D3  +  A2D2  +  A-jD  +  A0  =  0  (Al) 

Here  A^,  A .  .  .  Aq  are  functions  of  K,  £,  and  odq  for  a  flight 

condition  in  which  the  airplane  parameters  in  equations  (3)  to  (5)  are 
known.  If  any  one  of  the  three  autopilot  parameters  is  fixed,  stability 
boundaries  may  be  obtained  in  the  plane  defined  by  the  other  two  param¬ 
eters  by  letting  D  =  R  +  in},  fixing  the  value  of  R  for  each  curve, 
and  varying  cd.  Since  equation  (Al)  is  a  complex  equation,  it  may  be 
written  as  two  real  equations  and  solved  for  the  two  autopilot  variable 
parameters  at  each  value  of  D  =  R  +  in). 

According  to  reference  6,  if 


fn(R)=-i.^M  (n  =  1,2,3,'..  0  (A2) 

n*  dR 

then  the  two  real  equations  obtained  from  equation  (Al)  by  setting 
D  equal  to  R  +  icD  are 

F(R)  -  <n2f2(R)  +  uAf^(R)  -  cd6  =  0  (A3) 


fx(R)  -  o>2f3(R)  +  o>4f5(R)  =  0  (Ah) 

Since  the  coefficients  of  F(R)  and  fn(R)  are  functions  of  the 
two  variable  autopilot  parameters  only,  the  values  of  these  two  param¬ 
eters  which  will  yield  an  oscillatory  mode  of  motion  with  a  given 
damping  may  be  obtained  by  choosing  a  value  for  R  and  solving  equa¬ 
tions  (A3)  and  (a4)  at  any  value  of  cd.  This  procedure  yields  a  point 
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on  the  given  R  curve  in  the  parameter  plane;  for  instance,  in  the 
£ci>0-plane  if  £  and  cd0  he  chosen  as  the  variable  parameters.  In  this 
case,  calculations  are  simplified  by  the  use  in  equations  (A3)  and  (a4) 
of  A  and  B  as  defined  in  equation  (l),  since  these  equations  become 
linear  in  A  and  B . 


Because  of  the  increasing  difficulty  of  obtaining  accurate  values 
of  stability  derivatives  for  present-day  airplanes,  the  characteristics 
of  the  airplane  may  be  given  in  terms  of  transfer  functions  which  are 
obtained  from  flight  data.  In  order  to  obtain  the  form  of  the  char¬ 
acteristic  equation  in  terms  of  transfer  functions,  equation  (7)  is 

divided  through  by  (d2  +  2£<d0D  +  o>02)fo(d)  .  The  resulting  equation  is 


Kh>  2 

O 

V2  „  (&V>  -  l  %)  (^Kx^2  -  ^  c  !pD)  +  g  CLC 

D2  +  2  £<d0D  +  co02 

b2  Cn5  F„(D) 

The  first  factor  in  the  second  term  of  this  equation  is  the  autopilot 


transfer  function 


JP't  p 


(D)  (See  eq.  8.)  If  equations  (3)  to  (5)  are 


solved  operationally  for  zero  initial  conditions,  the  expression  in 
brackets  is  found  to  be  the  airplane  transfer  function: 


Dt 


JA 


(D)  = 


Vc 


"2  Un6 


SD  -  F  %)(2>W^  -  h  C*pD)  +  ^  CLCl, 

F0(D) 


(A5) 


Therefore,  the  characteristic  equation  may  be  written 


(D) 


(A  6) 


in  terms  of  the  autopilot  and  airplane  transfer  functions.  The  airplane 
transfer  function  for  the  lateral  motion,  as  shown  in  equation  (A5) , 
should  have  the  form  of  a  cubic  in  D  over  a  quartic  in  D.  Therefore, 
when  equation  (a6)  is  cleared  of  fractions  a  sixth-degree  equation  of 
the  form  of  equation  (Al)  is  again  obtained  as  the  characteristic 
equation. 
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The  setup  and  solution  of  equations  (A3)  and  (a4)  would  be  greatly 
simplified  if  the  characteristics  of  the  airplane  lateral  motion  could 
be  expressed  as  a  second-order  system  rather  than  a  fourth-order  system. 
Since  the  yaw  damper  is  used  essentially  for  damping  the  Dutch  roll 
oscillation,  it  was  thought  desirable  to  see  whether  the  airplane  could 
be  replaced  by  an  equivalent  oscillator  whose  damping  and  period  were 
those  of  the  Dutch  roll  oscillation  of  the  airplane.  That  is,  if  the 
Dutch  roll  mode,  as  calculated  when  three  degrees  of  freedom  are  con¬ 
sidered,  is  obtained  from  the  quadratic  equation  D2  +  PQD  +  Qq  =  0, 

the  question  is  whether  the  airplane  characteristics  may  be  represented 
by  the  oscillator  described  by 

-  gf)(°2  ♦  v  ♦  0  -  0 

Equation  (A5)  may  be  rewritten  in  the  form 


Dt 

8 


(D)  = 


V 

p?c%D 


K 


2MKz"  -  jp-)(D2  +  Po®  +  %) 


(D  +  si)(d£ 


+  a^D  +  a^ 


d(d2  +  a^D  + 


-C-jD 


D  +  PoD  +  Qo 


(D  +  al)( 


D  +  agD  + 


D  ^D2  +  aj|D  + 


(A7) 


'where  the  constant  C-^  is  defined  by 


V2C 


Ci¬ 


ne 


2b 2mJkz2  - 


K- 


X 


The  expression  in  the  first  bracket  on  the  right-hand  side  of  equa¬ 


tion  (A7)  is  the  same  as  would  be  obtained  for 


Dt 


lAja 


(D)  if  the  assump¬ 


tion  were  made  that  the  airplane  could  be  represented  by  the  equivalent 
oscillator  described  previously.  Thus,  if  the  expression  in  the  second 
bracket  has  a  frequency  response  which  does  not  significantly  differ 
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from  unity  for  the  range  of  frequencies  which  is  important  in  the  air¬ 
plane  frequency  response,  then  the  expression  in  the  first  bracket  can 
be  used  as  a  valid  approximation  to  the  airplane  frequency  response. 
Therefore,  the  expression  in  the  second  bracket  can  be  considered  as  a 
correction  function. 

The  factor  D2  +  a^D  +  a^  yields  the  two  real  characteristic  roots 

of  the  airplane  lateral  motion  corresponding  to  the  spiral  mode  and  the 
damping- in-roll  mode.  The  usual  approximation  to  the  damping- in-roll 

root  is  VC  Zp/^bii-bK^2 ,  and  the  spiral  root  is  of  the  order  of  yCyp^bp^. 

However,  the  numerator  of  the  second  bracket  can  be  shown  to  be 


v  fcb 

*\2Hb 


+ 


Czp  Vg  +  v£  °zpCyp 

H>Kx2/  *2  ®h»Kx2 


v3  CL% 

b3 


If  C2r  =  0,  this  cubic  has  the  roots 


D  =  0 


D  = 


J.  1 

2^  To 


C  -j' 

zp  V  , 

Vx2  b 


which  are  merely  approximations  to  the  real  roots  which  occur  in  the 
cubic  in  the  denominator  of  equation  (A7) •  It  can  therefore  be  expected 
that  the  second  bracket  will  usually  be  close  to  unity  and  that  no  large 
errors  will  arise  if  the  airplane  lateral  motion  is  represented  by  a 
single-degree-of -freedom  oscillator  in  yaw  with  nondimen sional  inertia 


2p-b|K22  -  J  and  the  same  period  and  damping  as  the  Dutch  roll  oscil- 

\  %  / 

lation.  This  is  especially  true  for  small  values  of  c2g. 


For  the  airplane  used  as  an  example  in  this  paper,  the  correction 
(D  +  a, ) (d2  +  a?D  +  a,) 

function  - - was  evaluated  throughout  the  fre- 

D^D2  +  a^D  +  a^y 

quency  range  and  its  effect  on  the  total  airplane  transfer  function  was 
found  to  be  very  small  except  close  to  zero  frequency.  Finally,  actual 
comparisons  of  the  stability  boundaries  for  the  three-degree-of-freedom 
analysis  and  the  equivalent-oscillator  approximation  showed  that  the 
approximation  was  valid. 
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When  equation  (A'j)  is  used  with  this  approximation,  the  character¬ 
istic  equation  as  given  in  equation  (A 6)  becomes 


1  + 


or 


-  S''^ 


K- 


■X 


K»o2 

C-jD 

D2,  +  2^0D  +  <nQ2 

D2  +  PqD  +  Qq 

+  PoD  +  «o)(d2  - 

2^qD  +  0>o2) 

=  0 


y2  2 

~  =  ° 


(A8) 


The  characteristic  equation  of  the  system  is  therefore  a  quartic, 

F(D)  =  D4  +  (PG  +  2£a>0)D3  +  (Qq  +  0>o2  +  2^(n0P0)D2  + 

^o^o  +  2^oQo  +  C  ]Kh>02^  D  +  Q0(n02 
=  D1*  +  (po  +  A)D3  +  (Qq  +  B  +  AP0)D2  +  1 
(poB  +  QqA  +  C^D  +  QqB 

=  0  (A9) 

The  equations  for  the  stability  boundaries  which  replace  equa¬ 
tion  (A3)  and  equation  (a4)  for  a  fourth-order  characteristic  equation 
are  , 

F(R)  -  o>2f2(R)  +  oS  =  0  (A10) 

fx(R)  -  cs2f3(R)  =0  ,  (All) 

The  functions  fn(R)  are  obtained  from  equation  (A2) . 
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As  an  example,  the  actual  parametric  equations  used  for  the  equiva¬ 
lent  oscillator  will  be  derived.  Equation  (A9)  is  used  in  equation  (A2) 
to  give 

fpCR)  =  ^R3  +  3R2(Pg  +  A)  +  2R(Qq  +  B  +  AP0)  +  (PQ  +  CpKjB  +  QqA 
f2(.R)  =  6r2  +  3R(P0  +  A)  +  Qq  +  PqA  +  B 
f3(R)  =  4R  +  PQ  +  A 


These  expressions  are  used  in  equations  (AlO)  and  (All) .  To  obtain  the 
constant-K  curves  (that  is,  curves  in  the  AB-plane  or  £aD0-plane) ,  the 
resulting  equations  are  written  in  A  and  B,  as  follows: 


3R  +  P0)<u2 


+  PqR  +  Qo 


+  (P0  +  CpK)R 


+  3PqR 


+  r2(r2 


+  PqR  + 


(A12) 


1 


(2R  +  PQ  +  CjK^B  = 


-(4R  +  Pq)cd2 


+  r(4Rc  +  3PqR 


(A13) 


For  a  given  value  of  K,  any  constant- damping  curve  is  obtained  by 
fixing  R  at  the  appropriate  value  and  taking  a  sufficient  number  of 
positive  values  for  co  to  obtain  the  points  necessary  to  determine  the 
curve.  The  resulting  two  linear  equations  in  A  and  B  are  solved  to 
obtain  each  point.  If  curves  in  the  £o>0-plane  are  desired  instead,  the 
values  of  £  and  tn0  can  be  obtained  for  each  set  of  values  of  A 
and  B  for  which  B  is  positive. 

The  constant-^  curves  in  the  Ko)0-plane  are  obtained  by  solving  the 
parametric  equations  for  K  and  tn0  instead  of  for  A  and  B.  For 
example,  solving  equation  (A13)  for  K  gives 
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K  =  j(A  +  p°  +  4r)“2  ■  [7po  +  2R)B  +  (Qo  +  2PqE  +  3R2)a 


r(2Qq  +  3PCR  +  4R2) 


+  po  +  ^R)^2  -  f(P0  +  2R)o>02  +  2£cd0(q0  +  PqR 

'1%  L  L 


k(2Qo  +  3P0R  +  4R2] 


+  3R2)  + 


(AHA) 


and  using  this  equation  in  equation  (A12)  gives 


(q0  -  R2  -  co2)  B  -  (PQ  +  2R)  Cos2  +.R2) A  +  o>^  -  ((^  +  2P0R  +  2E2)  (co2  +  R2)  = 

(qo  -  R2  -  (n2)o)02  -  2£cd0(p0  +  2r)  (co2  +  R2)  +  oA  -  +  2PQR  +  2R2)(cd2  +  R2)  =  0 

(A15) 

This  equation  is  a  quadratic  in  o>0  for  fixed  values  of  £  ,  R, 
and  0) .  The  values  of  o)0  from  equation  (A15)  are  then  used  in  equa¬ 
tion  (Al4)  to  obtain  the  points  on  the  curves  in  the  Ko^-plane.  Since 
equations  (A12)  and  (A13)  are  linear,  there  is  only  one  point  in  the 
AB-plane  or  the  £%-plane  for  each  oscillatory  mode  when  K  is  fixed. 

However,  for  a  fixed  value  of  £  there  may  be  two  points  in  the 

Kh>o -plane  for  a  given  oscillatory  mode  when  the  quadratic  equation  (A15) 

has  two  real,  positive  roots  for  u>0. 
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APPENDIX  B 

✓  )  *  ' 

DERIVATION  OF  EXPRESSIONS  FOR  MAXIMUM-DAMPING 
CHARACTERISTICS  WITH  FIXED  GEARING  RATIO 
'BY  USE  OF  EQUIVALENT -OSCILLATOR  CONCEPT 


The  characteristic  equation  of  the  airplane — yaw-damper  system, 
for  the  assumption  that  the  airplane  may  he  replaced  by  an  equivalent 
oscillator,  is 

D^  +  (po  +  a)d3  +  (Qq  +  B  +  AP0)p2  +  (PqB  +  AQ0  +  C^D  +  QqB  =  0 

(Bla) 

The  conditions  on  A  and  B  which  result  in  the  maximum  damping  for 
a  given  gearing  ratio  can  be  calculated  as  follows.  -For  the  case  of 
maximum  damping,  equation  (Bla)  will  have  two  pairs  of  equal  complex 
roots;,  hence  this  equation  becomes: 

(d2  +  PD  +  q)2  =  D^  +  2PD3  +  (p2  +  2q)d2  +  2PQD  +  Q2  (Bib) 

2 

where  D  +  PD  +  Q  is  the  quadratic  corresponding  to  the  double  oscil¬ 
latory  root  at  the  cusp  point.  The  following  four  equations  are  obtained 
by  equating  like  coefficients  in  equation  (Bib): 


2P  =  PQ  +  A 

(B2) 

P2  +  2Q  =  Qq  +  B  +  AP0 

(B3) 

2PQ  =  P0B  +  AQq  +  CjKB 

(B4) 

Q2  =  QqB  ' 

(B5) 

Since  K,  C]_,  PQ,  and  Q0  are  known,  these  four  independent  equations 

pay  be  solved  for  A,  B,  P,  and  Q.  Thus,  the  autopilot  parameters 
which  will  give  the  highest  damping  to  the  Dutch  roll  oscillation,  and 
also  the  period  and  damping  of  this  oscillation,  may  be  determined. 
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From  equations  (B2)  and  (B5) 

A  =  2P  -  P0 


(b6) 


(B7) 


Using  equations  (b6)  and  (B7)  in  equation  (B4)  yields 


(P0  -h  CxK)Sg  -  P0Q0g 
2«o(«  -  «o) 


(E8) 


When  equations  (b6)  and  (B7)  are  used  in  equation  (B3) ,  the  result  may¬ 
be  written  - 


«o(p  -  P°f  -  («  -  «o f 


(B  9) 


and  by  substituting  equation  (b8)  into  equation  (B9)  a  quartic  equation 
in  Q  is  obtained,  which  may  be  written 


Q0)  +  C-lKQ' 


H>(Q  -  %)k 


Finally,  this  equation  yields  the  two  quadratics  in  Q: 


- )  Q2  -  2Q0Q  +  Qq2  =  0  (BIO) 

Po  ±  2^/ 

Evaluation  of  the  discriminant  of  equation  (BIO)  reveals  that  equa¬ 
tion  (BIO)  has  only  two  real  roots,  one  less  than.  Q0  and  the  other 

greater  than  Qq.  As  can  be  seen  from  figure  6,  the  smaller  value  of  Q 

(which  yields  a  smaller  value  of  o>0)  corresponds  to  the  breakdown  of  the 
unstable  loops  for  positive  K.  The  larger  real  root  of  equation  (BIO) 
is  therefore  the  one  which  yields  the  maximum-damping  point  for  positive 
gearing  ratio.  Negative  gearing  ratios  will  be  discussed  subsequently. 

The  corresponding  value  of  P  is  obtained  by  substituting  this 
larger  root  into  equation  (B8) ,  and  the  maximum- damping  point  in  the 
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AB-plane  is  easily  obtained  from  equations  (B6)  and  (B7) •  Thus,  the 
value  of  the  maximum  damping  of  the  Dutch  roll  oscillation  obtainable 
with  a  given  gearing  ratio,  'and  also  the  values  of  autopilot  frequency 
and  damping  which  will  yield  this  maximum  damping,  may  be  obtained  by 
simply  solving  a  quadratic  equation  (eq.  (BIO)). 

The  procedure  is  even  simpler  when  the  required  damping  is  given 
and  the  values  for  the  autopilot  parameters  which  would  make  this 
damping  the  maximum  damping  are  to  be  determined.  The  value  of  P  is 
determined  by  the  required  damping,  since 


1-386 

Tl/2 


(Bll) 


and  the  corresponding  optimum  value  of  A  may  be  obtained  from  equa¬ 
tion  (B6) .  Equation  (B9)  now  gives 


Q  =  Q0  ±  ^  (P  -  P0) 


(B12) 


Both  values  of  Q  obtained  from  equation  (B12)  correspond  to  optimum 
points,  since  they  both  correspond  to  the  given  positive  value  of 
damping.  These  values  may  now  be  used  in  equation  (B7)  to  obtain  the 
corresponding  values  of  B.  Finally,  the  two  values  of  gearing  ratio 
may  be  obtained  from  equation  (B8) ,  which  gives 


2PQ  -  P0B  -  AQ0 


(B13) 


The  smaller  values  of  ,  Q  (obtained  from  using  the  minus  sign  in 
eq.  (B12) )  result  in  negative  values  of  gearing  ratio  and  values  of  o>0 
less  than  the  airplane  frequency.  The  possibility  of  using  negative 
gearing  to  improve  the  damping  is  discussed  further  in  the  body  of  the 
paper . 
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TABLE  I 

PARAMETERS  OF  AIRPLANE  USED  AS  EXAMPLE  TO  ILLUSTRATE 
USE  OF  CONSTANT -DAMPING  CURVES 


Altitude,  ft  ....  . 
Wing  loading,  lb/sq  ft 

V,  ft/sec  . 

b,  ft . . 

CL  .........  . 

^b  •  * . * 


KXZ . 

Cjp?  Per  radian 

C 2r,  per  radian 
Cnp>  Per  radian 
Cn;r,  per  radian 
Cy^,  per  radian 
Cnp,  per  radian 
Cjp,  per  radian 
Cng,  per  radian 
Period,  sec  .  .  . 
Tl/2'  sec  •  •  •  • 

PQ,  sec-"*-  :  .  .  . 

-2 

Q0,  sec  •  •  •  . 

_  -2 

Cp,  sec  .  .  .  . 


30,000 
.  .  65 

.  .  797 
.  .  28 
.  0.23 

.  80.7 

0.0097 

0.0513 

-0.00145 

.  -0 . 4o 

0.80 
.  -0.02 
.  -0.40 

-1.0 
0.25 
.  -0.126 
.  -0.163 

.  .  1.3 

.  .  2.6 

.  0.537 

.  23.84 

.  15.98 
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TABLE  II 

PARAMETERS  USED  FOR  THREE  TYPICAL  FLIGHT  CONDITIONS  OF  THE 
AIRPLANE  IN  DESIGN  OF  COMPROMISE  YAW-RATE  DAMPER 


Altitude,  ft  .  . 

Case  I 

Sea  level 

Case  II 

30,000 

Case  III 

30,000 

Mach  number  . .  . 

. 

• 

• 

0.22 

0.75 

0.5 

Wing  loading,  lb/sq 

ft 

5b 

65 

85 

^b . .  * 

• 

• 

25.2 

81.25 

106.3 

Kx2  ...... 

• 

• 

• 

• 

• 

0.008l 

O.OO69 

0.0051 

KZ2 . 

• 

• 

• 

• 

• 

• 

• 

0.0433 

0.0419 

0.0409 

Kxz . . 

• 

■ 

• 

• 

• 

• 

• 

0.0027 

0.0025 

0.0010 

CL  ......  . 

• 

• 

• 

• 

• 

• 

• 

■ 

0.765 

0.262 

0.771 

Cnp . 

• 

• 

• 

• 

• 

• 

• 

0.205 

0.205 

0.212 

% . 

• 

• 

• 

• 

• 

• 

• 

-0.099 

-0.107 

-0.095 

cYp  . . 

• 

• 

• 

• 

• 

• 

• 

-0.930 

-0.878 

-0.884 

Cz  . 

• 

• 

-0.425 

-0 . 474 

-0.435 

C2  . 

• 

■ 

• 

0.288 

0.200 

0.300 

cn  ...... 

np 

• 

% 

• 

0.003 

0.010 

0.003 

0nr 

• 

• 

• 

-0 . 165 

-0.150 

-0.165' 

% . 

• 

• 

• 

-0.163 

-O.163 

-O.163 

Period,  sec  .  . 

• 

2.3 

1.4 

2.1 

■T1/2,  sec  ... 

• 

2.0  ' 

6.9 

2.4 

PQ,  sec  ... 

« 

0.704 

0.200 

,  0.573 

Qq,  sec'2  .  .  . 

• 

7-79 

21.4 

8.78 

C-^,  sec-2  .  .  . 

• 

5.63 

17.0 

5.92 

v/b,  sec"-'-  .  .  . 

• 

• 

. 

8.68 

26.64 

17.77 

second  are  given  at  representative  points 
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Figure  8.-  Comparison  of  motions  with  three  different  autopilots  whose  characteristics 
on  a  curve  for  Tw2  =  0.60  second  (corresponding  to  points  1,  2,  and  3  in  fig.  6). 
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Figure  11.-  Values  of  autopilot  gearing  ratio  necessary  to  obtain  a 
given  damping  as  the  maximum  damping. 


radians  /sec 
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16.-  Values  of  autopilot  parameters  K  and  odq  for  which  an 
llation  will  occur  in  the  lateral  motion  with  ~  0.60  second. 

0.3.  Values  of  to  in  radians  per  second  are  given  at -representative 
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20.-  Comparison  of  the  boundary  of  equal  roots,  the  airplane- 
ing  region,  and  a  typical  increased -damping  region  in  the 
plane  for  values  of  £  *  0.3>  0.6,  and  0.9- 


subsequent  to  a  5°  sideslip  disturbance.  Autopilot 
K  =  0.14  deg/deg/sec:  £  =  0.523}  u>o  =  9-49  radians  ; 


ure  24.-  Calculated  motions  for  three  flight  conditions  of  the  airplane 
subsequent  to  a  5°  sideslip  disturbance.  Autopilot  characteristics: 

K  =  0.12  deg/deg/sec;  £  =  0.485;  =  8.8l  radians  per  second. 
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